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1. Introduction. This paper can be seen as part of an ongoing project 
whose aim is to give a pathwise definition to stochastic PDEs. Indeed, the 
rough path theory [5, 13, 16, 17] and its variants [4, 6] have now reached a 
certain level of maturity, leading to a proper definition of differential equa- 
tions driven by irregular signals and in particular by a fractional Brownian 
motion [2]. Starting from this observation, we have tried in [12] to define 
and solve the following general problem: let B be a separable Banach space, 
and A:D(A) —¥ B the infinitesimal generator of an analytical semigroup 
{S t ; t > 0} on B, inducing the family {B a ;a G R} with B a = D{(-A) a ). Let 
also / be a function from B to C(B- a ,B~ a ) for a given a > and x a noisy 
input, considered as a function from M + to B- a . Then, for T > 0, consider 
the equation 



with an initial condition y$ G B. The main example we have in mind is 
the case of the 1-dimensional heat equation in [0,1], namely B = L 2 ([0, 1]), 
A = A with Dirichlet boundary conditions, the usual Sobolev spaces B a = 
H a = W 2a > 2 , and x a fractional Brownian motion with Hurst parameter H 
taking values in B- a . Notice, in particular, that we wish to consider a noise 
x which is irregular in both time and space. Then, in [12], we gave a local 
existence and uniqueness result for equation (1), by considering it in its mild 
form 



where we let Sts = St- S and interpreting the integral in this mild formula- 
tion as a Young integral. Once the equation is set under the form (2), the 
main problem one is faced with is to quantify the regularization of the semi- 
group Sts on the term f(y s ) dx s , and then to elaborate the right fixed point 
argument in order to solve the equation. The general results of [12] could 
be applied in the case of the stochastic heat equation driven by a fractional 
Brownian motion with Hurst parameter H > 1/2. They should be compared 
with the reference [18], where a nonlinear fractional SPDE is solved thanks 
to some fractional calculus methods, but where x is a smooth noise in space. 

In the current article, we would like to go one step further with respect 
to [12], and set the basis of a real rough path expansion in order to define 



(1) 



dy t = Ay t dt + f(y t ) dx t 



te[0,T\ 
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and solve equation (2), which would allow to consider, in the case of the 
heat equation in [0,1], a fractional Brownian motion with Hurst parameter 
H < 1/2. This task is quite long and involved, but let us summarize at this 
point some of the ideas we have followed. 

(1) We will recast equation (2) in a suitable way for expansions accord- 
ing to the following simple observation: we have tried to solve our evolution 
equation by means of its infinite-dimensional setting, since it allows to con- 
sider x and y as functions of a unique parameter t G [0, T], which makes 
its rough path type analysis easier (see [11] and [22] for a multiparametric 
setting). However, when we come to the applications to the heat equation, 
we will consider the evolution equation in [0,T] x [0, 1] under the form 



where G stands for the fundamental solution to the heat equation, a : R — >■ 
R is a regular function, and x{ds,drj) is understood as the distributional 
derivative of a real- valued continuous process on [0, T] x [0, 1] . This definition 
of our equation is of course equivalent to (2) when / is considered as the 
pointwise nonlinear operator [f(yt)](0 = ^{UtiO)- Now, when written under 
its multiparametric form (3), the equation is also equivalent to 



and it happens that this simple reformulation is much more convenient for 
our future expansions than the original one. When we go back to the original 
infinite-dimensional setting, we can recast (2) into 



where / is now a smooth function from B to B, and x will be understood 
as a Holder-continuous process taking values in a space of deregularizing 
operators from B to a distributional space B_^ for a certain Q > 0. The 
product dx s f{y s ) will then be regularized again by the action of St s , in 
a way which will be quantified later on. Notice that the form (4) of our 
evolution equation is a little unusual in the SPDE theory, but makes sense 
in our context. 

(2) Instead of considering Riemann sums like in [12] or like in the original 
Lyons' theory [16], our analysis will be based on the theory of generalized 
differentials, called /c-increments, contained in [6]. Roughly speaking, this 
theory is based on the fact that an elementary operator, called 5, can trans- 
form an integral f l dg u [h u — h s ], seen as a function of the variables s and t, 
into a finite difference product (gt — g s ){ht — h s ). Furthermore, under some 
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additional regularity properties on g and h, the operator 5 can be inverted, 
and its inverse A, called sewing map (from [4]), will be the building stone 
of our extension of the notion of integral. Notice that, whenever g and h 
are Holder-continuous with Holder exponent > 1/2, this extension coincides 
with the usual Young integral. When we consider an integral of the form 
J* dg u 4>(g u ) for a Holder-continuous function g with Holder exponent in 
(1/3,1/2] admitting a Levy area, our definition of integral also coincides 
with Lyons' one, as shown in [6]. In fact, if the usual rough path theory 
gives a richer point of view on the algebraic structure of the path x, it is 
worth mentioning that our approach has at least two advantages: 

(1) Once our unusual setting is assimilated, it becomes quite easy to figure 
out how a given expansion in terms of x can be leaded. And indeed, it 
will become clear throughout the paper, that the /c-increments theory 
provides a tool allowing some natural computations for our generalized 
integrals. 

(2) The only step where a discretization procedure is needed is the con- 
struction of the A map alluded to above, and this avoids some of the 
cumbersome calculations which are one of the main ingredients of the 
rough path theory. 

We hope that this paper will advocate for the use of the ^-increments theory, 
which obviously does not exclude the other approaches [4, 16]. 

(3) The fact that we are dealing with an evolution problem will force us 
to change some of the algebraic structure we will rely on, especially if one 
wants to take advantage of the regularizing effect of St- This will lead us to 
introduce an operator a ts = St s — Id for t > s, and a modified 5 operator, 
called 5, defined by 5 = 5 — a. The whole increment theory will have to be 
built again based on this modified operator, and we will see that it is really 
suitable for the evolution setting induced by (4). In particular, we will be able 
to define analogs of the Levy area and of the higher-order iterated integrals, 
which are of course harder to express than in the finite-dimensional case, 
but can be written, in the bilinear case [that is cr(r) = r in (3)], as 



Obviously, a convenient definition of iterated integrals is the key to reach 
the case of a Holder-continuous noise of order < 1/2. 

(4) The whole integration theory can be expressed in an abstract way, by 
just supposing a certain set of assumptions on some incremental operators 
like A 2 and A 3 . However, we will try to check these assumptions in some 
interesting cases, like the infinite-dimensional fractional Brownian motion for 
our Young type integration, or the infinite-dimensional Brownian motion for 




etc. 
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our step 3 expansion, based on X 2 and X 3 . Notice that the rough expansions 
for the fractional Brownian motion should be investigated in details too, 
but one is faced with an additional problem in this situation: on one hand, 
a Stratonovich type integration requires a lot of regularity in space for the 
noise, due to the well-known presence of some trace terms. On the other 
hand, the Skorokhod integral does not fulfill the algebraic requirements we 
ask for our integral extension. A discussion of these problems and some 
ideas to solve them will be included at the end of the paper, but for sake 
of conciseness, we will postpone a complete development of this part to a 
subsequent paper, and stick here to the Brownian case. 

This paper is structured as follows. In Section 2, we recall the basic setup 
of [6] which allows to embed the theory of rough paths in a theory of integra- 
tion of generalized differentials, called here k-increments. We wrote it with 
the aim of having a self-contained and pedagogical introduction to the topic. 
However, we give also a new and very elementary proof of the existence of 
the basic integration map A of [6]. In Section 3, we introduce and study a 
modified coboundary induced by the operator 5 on the complex of incre- 
ments, using the additional data provided by an analytic semigroup S, in 
such a way that the new complex can be shown to act simply on convolution 
integrals of the form appearing in equation (4) and on their iterated versions. 
This new complex maintain many of the properties of the original complex 
(e.g., its cohomology is trivial) and it is shown that when equipped with 
Holder-like norms which measures "smallness" of the increments, it admits 
a map, called A here, which is the main tool for building an integration (or 
better, convolution) theory over those 1-increments which are good enough 
(again, in a suitable sense, to be specified in due time). A key feature of this 
perturbed complex is that, due to the convolution with the semigroup S, 
"space" and "time" regularity of increments depends on each other: we can 
gain space regularity by loosing some time regularity and vice-versa. This 
property will be essential for the solution of the evolution problem by fixed- 
point arguments. In Section 4.2, we use the theory outlined in Section 3.2 to 
define the convolution integral in the Young sense and solve a class of non- 
linear evolution problems, reobtaining some results of the work [12]. Notice 
that we will also improve some of our previous results contained in [12], in 
the sense that we will be able to construct global solutions to our evolution 
equations in the Young context. In Section 5, we study the bilinear evolution 
problem 



(6) 




We will also introduce a notion of rough-path suitable for noises driving 
evolution equations. By exploiting this pathwise technique, we are able to 
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obtain automatically the flow semigroup of the equation and we will show 
how to express this semigroup as a convergent series of iterated-integrals 
which are the lift of the step-3 rough path used in the construction of the 
solution. In Section 6, we turn to a nonlinear case of evolution system, 
namely the case of the quadratic type equation 



where B stands for the pointwise multiplication of functions. This requires 
the additional careful introduction of a collection of a priori increments 
indexed by planar trees, and an associate notion of controlled path. Finally, 
all our results will be applied in the concrete case of the stochastic heat 
equation on the circle, in a setting recalled at Section 3.4. The case of a 
fractional Brownian case is handled in the special situation of the Young 
theory, while we stick to the example of an infinite-dimensional Brownian 
motion in the rougher situation. We build the rough path associated to this 
latter noise and provide concrete conditions where the theory outlined in 
the previous sections can be fruitfully applied. A systematic study of the 
regularity properties of the incremental operators defined as X 2 or X 3 in 
(5) will also be provided at Sections 6.5 and 6.6, thanks to some Feynman 
diagram techniques. 

2. Algebraic integration in one dimension. The integration theory intro- 
duced in [6] is based on an algebraic structure, which turns out to be useful 
for computational purposes, but has also its own interest. Since this setting 
is quite nonstandard, compared with the one developed in [16], and since 
we will elaborate on it throughout the paper, we will recall briefly here its 
main features. We also provide an elementary proof of the existence of the 
A map. 

2.1. Increments. As mentioned in the Introduction, the extended inte- 
gral we deal with is based on the notion of increment, together with an 
elementary operator d acting on them. However, this simple structure gives 
rise to a nice topological structure that we will describe briefly here: first 
of all, for an arbitrary real number T > 0, a vector space V, and an inte- 
ger k > 1, we denote by Cfe(V) the set of functions g: [0,T] k — > V such that 
gtx—u = whenever t{ = tj+i for some i < k — 1. Such a function will be 
called a (k — 1) -increment, and we will set C*(V) = Ufc>i^fc(^)- The oper- 
ator S alluded to above can be seen as a coboundary operator acting on 
/c-increments, inducing a cochain complex (C*,<5), and is defined as follows 




on C k (V): 



k+l 



(7) 



8:C k {V)^C k+1 {V) 



(fa)ti-t, 



: fe+i 



i=l 
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where tj means that this particular argument is omitted. Then a funda- 
mental property of 5, which is easily verified, is that 55 = 0, where 55 is 
considered as an operator from Cfc(V) to Ck + 2(V). We will denote ZC^iV) = 
C k (V) n Ker(5| Cfc(y) and #C fc (V) := C k (V) n Im5| Cfc respectively the 

spaces of k-cocycles and of k-coboundaries, following standard conventions 
of homological algebra. 

Some simple examples of actions of 5, which will be the ones we will really 
use throughout the paper, are obtained by letting g G Ci and h G 62- Then, 
for any t,u,s G [0, T] , we have 

(8) {5g)ts=9t-g s and (5h) tus = h ts - h tu - h us . 

Furthermore, it is readily checked that the complex (C*,<5) is acyclic, that 
is, ZCk+iiV) = BCk{V) for any k > 1, or otherwise stated, the sequence 



is exact. In particular, the following basic property, which we label for further 
use, holds true. 

Lemma 2.1. Let k>l and h G ZCk+i(V) . Then there exists a (nonunique) 
f G Cfc(V) such that h = 5f . 

Proof. This elementary proof is included in [6]; see also Proposition 3.1 
below. Let us just mention that ft v -t h = ht v ..t k o 1S a possible choice. □ 

Remark 2.2. Observe that Lemma 2.1 implies that all the elements 
h G C2(V), such that 5h = 0, can be written as h = 5f for some (nonunique) 
/ G C\(y). Thus, we get a heuristic interpretation of <5|c 2 (v) : it measures how 
much a given 1-increment is far from being an exact increment of a function 
(i.e., a finite difference). 

Notice that our future discussions will mainly rely on fc-increments with 
k < 2, for which we will use some analytical assumptions. Namely, we mea- 
sure the size of these increments by Holder norms defined in the following 



0) 



-> R -> Ci(V) A C 2 {V) A C 3 (V) A C 4 (V) -»••■• 



way: for / G C 2 (V) let 




and C£(V) = {f€C 2 (V);\\f\\ lt <oo}. 



In the same way, for h G Cs(V), set 



INI 



sup 

s,u,te[0,T] 



u - spit - u\p 



(10) 
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where the last infimum is taken over all sequences {hi £ C3(V)} such that 
h = Yli^i and for all choices of the numbers pi G (0,/u). Then || • ||^ is easily 
seen to be a norm on C^iV), and we set 

c 3 M 0/) := {/i e c 3 (niNU < °o}. 

Eventually, let C\ + {V) = Ujt>l C 3 00> and remark that the same kind of 
norms can be considered on the spaces ZCziV), leading to the definition of 
some spaces ZC%{V) and ZC\ + {V). 

With this notation in mind, the following proposition is a basic result 
which is at the core of our approach to pathwise integration. 

Proposition 2.3 (The sewing map A). There exists a unique linear 
map h: ZC\ + (V) -^-C\ + (y) (the sewing map) such that 

6A = Id 2 c 3 (V)- 

Furthermore, for any fi>l, this map is continuous from ZC^(V) to C^V) 
and we have 

(11) im^^z^WH^ heZC l 3 + (V). 

Proof. For sake of completeness, we include a proof of this result here, 
which is more elementary than the one provided in [6], and which will be 
generalized at Theorem 3.5. For the sake of notation, we will omit the de- 
pendence in V in our functional spaces, and write for instance C3 instead of 
C^(V). Let then h be an element of ZC% C ZC\ + for some fj, > 1. 

Step 1 : Let us first prove the uniqueness of the 1-increment MsC^ 1 such 
that SM = h. Indeed, let M, M be two elements of C% satisfying SM = 5M = 
h and set Q = M — M. Then 5Q = and Q € • Invoking Lemma 2.1, there 
exists an element q 6 C\ such that Q = 5q, but since fj, > 1, q is a function 
on [0,T] with zero derivative, that is a constant and then Q = 0. 

Step 2: Let us construct now a process MgCj, with fi>l, satisfying 
SM = h. Since 5h = 0, invoking again Lemma 2.1, we know that there exists 
a B G C2 such that 8B = h. Pick s,t G [0,T], such that s < t in order to fix 
ideas, and for n > 0, consider the dyadic partition {r™; i < 2 n } of the interval 
[s,t], where 

(12) = s + for i < 2 n . 
Then for n > set 

2 n -l 

(13) M£ = Bu-Y, B ^rt- 

i=0 
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Then it is readily checked that M® s = 0. Furthermore, we have 

2 n -l 

MJ^ 1 — M™ = /J (B n+1 n+1 — B n+1 n+1 — £ n+1 n+1 ) 
— ' r 2i+2> r 2i 7 2i+l' r 2i r 2i+2> T 2i+l 

j=0 

2 n -l 2™-l 

= / . (SB) n + l n + 1 n+1 = > /l„n + l n + 1 n+1, 
^— ' '2i+2>' 2i+l>'2< z — ' ' 2i+2' r 2i+l''2i 

j=0 i=0 

and since /i £ C3 with /x > 1 , we obtain 

i M n _ )M-n+l| < INUfo ~ g ) M 
A2 ta I S 2"0*-l) ' 

which yields that Mt s = lim n _ ) . 0O M™ s exists, and satisfies inequality (11). 
Step 3: Let us consider now a general sequence {7r n ;n > 1} of partitions 

K > *i , • • • » r L > of t s > *] » with s = r o < r i < • • • < r t < r k n +i = *• We 

assume that 7r ra C vr ra+ i, and linij^oo k n = 00. Set 

(14) K B =£fc-£% + i.T- 

i=o 

It is easily seen that there exists 1 < / < k n such that 

2\t-s\ 



(15) K +l -rU\< , • 

Pick now such an index I, and let us transform 7r n into 7r, where 

Then, as in the previous step, 



Ml = Ml" - {5B) rT+iKirU = M*« - hr^n^, 



using the definition of the space C3 and the bound (15) we have 



|A/£-Af£"| <2** 



Repeating now this operation until we end up with the trivial partition 
7To = l° r which Mjj° = 0, we obtain 

k n OO 

\mz?\ < 2»\\hut - srj^r" < ^\\h\i\t - s\»Y,^ = c ^\t - sr. 

3=1 j=i 

Hence, there exists a subsequence {7r m ; m > 1} of {n n ;n > 1} such that M^ m 
converges to an element Mt s , satisfying Mt s < Cnh\t ~ S I M - With the same 
considerations as in [13], it can also be checked that the limit M does not 
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depend on the particular sequence of partitions we have chosen, and thus 
coincides with the one constructed at Step 2. 

Step 4 '■ It remains to show that SM = h. Consider then < s < u < t <T, 
and two sequences of partitions ir™ s and tt? u of [s,u] and [u,t], respectively, 
whose meshes tend to as n — >■ oo. Set also 7r" s = 7r^ U7r" s . From the previous 
step, one can construct easily some subsequences vr^,7r™,7r™, with 7r™ = 
<UtC, such that 

■7T m m TT m 

lim M^ u = M tu , lim Ml r = M us , hm = M ts . 

m— too m— ¥oo m— >oo 

Call now k™ (resp. kp^,k^ s ) the number of points of the partition 7r£? (resp. 
7r^,7r™). Then a direct computation, using definition (14), shows that for 
any < i < 2 n we have 




1 n-i' i / ^ ' i ' i+i 



Taking the limit m — > oo in the latter relation, we get (5M)t us = htusi which 
ends the proof. □ 



We can now give an algorithm for a canonical decomposition of the preim- 
age of the space ZC\ + (V), or in other words, of a function g G C2OO whose 
increment 5g is smooth enough. 

Corollary 2.4. Take an element g G C2CV), such that Sg G C3 (V) for 
/_i > 1. T/ien 5 can be decomposed in a unique way as 

g = 5f + A5g, 

where f G C\(V). 



PROOF. Elementary; see [6]. □ 



At this point, the connection of the structure we introduced with the 
problem of integration of irregular functions can be still quite obscure to the 
noninitiated reader. However, something interesting is already going on and 
the previous corollary has a very nice consequence which is the subject of 
the following property. 
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Corollary 2.5 (Integration of small increments). For any 1-increment 
g G C2(V), such that 8g € C\ + , set 5f = (Id — A5)g. Then 

n 

{Sf) ts = Jim y2g u+lU , 

Jits — M) . 

where the limit is over any partition lLj S = {^o = ^> • • • > = s } of [t, s] whose 
mesh tends to zero. The 1-increment 5f is the indefinite integral of the 1- 
increment g. 

Proof. Just consider the equation g = 5f + A5g and write 

n n n 

i=0 i=0 i=0 

n 

= (*/)*+ £(A^W 
i=0 

Then observe that, due to the fact that A5g 6 C2 + (l^), the last sum converges 
to zero. □ 

2.2. Computations in C*. For sake of simplicity, let us assume, until Sec- 
tion 3, that V = R, and set C&(R) = Ck- Then the complex (C*, 5) is an (as- 
sociative, noncommutative) graded algebra once endowed with the following 
product: for g GC„ and /i £ C m let gh € C n+m _i the element defined by 

(16) {gh)t- i .,...,t m+n -i =3ti,...,tn^tn,...,Wn-l! *1) • • • iWn+1 6 [0,T]. 

In this context, the coboundary 5 act as a graded derivation with respect to 
the algebra structure. In particular, we have the following useful properties. 

Proposition 2.6. The following differentiation rules hold true: 

(1) Let g,h be two elements of C\. Then 

(17) 5(gh)=5gh + g5h. 

(2) Let g € C\ and h^Ci- Then 

S(gh) = 5g h + g Sh, 5(hg) = 8hg — h5g. 

Proof. We will just prove (17), the other relations being equally trivial: 
if g, h G Ci, then 

[5(gh)] ts =gth t -g s h s = gt{h - h s ) + (g t - g s )h s = g t (5h) ts + (6g) ts h s , 
which proves our claim. □ 
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The iterated integrals of smooth functions on [0,T] are obviously partic- 
ular cases of elements of C which will be of interest for us, and let us recall 
some basic rules for these objects: consider f,g G Cf 3 , where Cf 3 is the set 
of smooth functions from [0,T] to R. Then the integral J dgf, which will 
be denoted by J[dg /), can be considered as an element of C|°. That is, for 
s,t € [0,T], we set 

t 



Jts{dgf)= (^J dgf^j = J dg u f u . 



The multiple integrals can also be defined in the following way: given a 
smooth element h £ C|° and s,t 6 [0, T], we set 



J ts {dgh)= (^J dgh^j = J dg u h us . 



ts 

In particular, the double integral Jt s {df^ df 2 f 1 ) is defined, for f ,f 2 ,f 3 € 
, as 

Jts(dfdf 2 f 1 )=(J dfdff^J = ^df'Jusidff 1 ) 
and if f\...,f n+1 eCf, we set 

(18) Jts{df n+l dr---df 2 f l ) = j\f n a +l J us {df n ---df 2 f 1 ), 

which defines the iterated integrals of smooth functions recursively. 

The following relations between multiple integrals and the operator 5 will 
also be useful in the remainder of the paper. 

Proposition 2.7. Let f,g be two elements o/Cf°. Then, recalling the 
convention (16), it holds that 

5f = J(df), 5(J(dgf)) = 0, 

5(J(dgdf)) = (6g)(5f) = J(dg)J(df) 

and, in general, 

rt-l 

5(J(df n ■ ■ ■ df 1 )) = Y,J( d f n ■ ■ ■ df +1 )J(df ■ ■ ■ df 1 ). 

i=l 

Proof. Here again, the proof is elementary, and we will just show the 
third of these relations: we have, for s,t G [0, T], 

Jts (dg df) = / dg u (f u - f s ) = / dg u f u - K ts , 

J S J S 
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with Kts = (gt — 9s) fs- The first term of the right-hand side is easily seen to 
be in ZCi- Thus, 

S(J(dgdf)) tus = -(5K) tus = [g t - g u ][f u - f s ], 

which gives the announced result. □ 

2.3. Dissection of an integral. The purpose of this section is not to pro- 
vide an account on all the computations contained in [6]. However, we will 
go into some semi- heuristic considerations that, hopefully, will shed some 
light on the way we will solve rough PDEs later on: with the notation of 
Section 2.2 in mind, we will try to give, intuitively speaking, a meaning to 
the integral j (p(x) dx = J[dx <p(x)) for a nonsmooth function x £C\. Notice 
that, in the sequel, x should be considered as a vector valued function, since 
the whole theory can be handled via the Doss-Soussman methodology in the 
real case. However, we will present the main ideas of the algorithm below as 
if x were real valued, the generalization from K to M n being just a matter 
of (cumbersome) notation. 

2.3.1. The Young case. The first idea one can have in mind in order 
to define J~(dx(p(x)) is to perform an expansion around the increment dx: 
indeed, in the smooth case, we have 

(19) J(dx tp(x)) = 5x <p(x) + J(dx dip(x)). 

If we wish to extend the right-hand side of (19) to a nonsmooth case, we 
see that the first term is harmless, since it is defined independently of the 
regularity of x, by 

[5x(p(x)] ts = [x t - x s ]ip(x s ) for s,te [0,T]. 

The last term of (19) is more problematic and we proceed to its dissection 
by the application of 5: invoking Proposition 2.7, we get, in the smooth case, 
that 

5{J(dxdif{x))) = 8x8((p(x)), that is, 

(20) 

[d(J(dxd(p(x)))]tus = [Sx] tu [S(ip(x))] us . 

Now the r.h.s. of (20) is well defined independently of the regularity of 
x. Thus, if 5x5(f(x)) 6 C\ + , which happens when x € Cf with a > \ and 
ip G C 1 (M), then Proposition 2.3 can be applied, and A[5x 5(tp(x))] is defined 
unambiguously. Hence, owing to (20), we set 

J(dxd<p(x))=A{SxS(<p(x))) 

and 



(21) J(dx(p(x))=6xip(x)+A(5x5((p(x))) = (Id - AS)[Sxip(x)], 
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where the last equality is due to Proposition 2.6 and to the fact that 55x = 0. 
Notice once again that this construction is valid whenever x 6 Cf with ol>\ 
and ip £ C 1 (1R), and it is easily shown, along the same lines as in the proof 
of Proposition 2.3 that the integral J{dxip{x)) defined by (21) corresponds 
to the usual Young integral. 

2.3.2. Case of a a-Holder path with ^ < a < ^. The construction (21) 

does not work if x ^ C\^ 2+ . However, if x € Cf with a > g, we can proceed 
further in the expansion of equation (19) by observing that, still in the 
smooth case, we have, for s,t £ [0,T], 

t pt ft fU 

[dip(x)] u = J dx u ip'(x u ) = [x t -x s ]ip'(x s )+ / dx u / dx v (p"(x v ), 

Js Js Js 

or according to the notation of Section 2.2, 

(22) 5ip(x) = J(d<p{x)) = J(dxip'(x)) = 5xip'(x) + J(dxdip'{x)). 
Injecting this equality in equation (19), thanks to (18), we obtain 

(23) J(dxip{x)) = 5xtp(x) + J \dx dx)tp' '(x) + J \dx dx dip' \x)) . 

Let us assume now that we are given a process J{dxdx) £ C2, usually (and 
somewhat improperly) called the Levy area of x, such that 

(24) 5(J(dxdx))=6x5x and J(dxdx)eC% a . 

This assumption is of course not automatically satisfied, but it can be 
checked for instance in the Brownian and fractional Brownian cases. Then 
the right-hand side of (23) is again well defined independently of the regu- 
larity of x, except for the last term. However, recast equation (23) as 

— J (dx dx dip' '(x)) = —3{dxip(x)) + 5xip(x) + 3{dx dx)ip' (x), 

and apply again 5 to both sides of this last expression. Invoking Proposition 
2.7 and recalling that 5(3 [dx dx)) = 5x5x, we obtain 

— 53 (dx dx dip' (x)) = —5x 5p(x) + 5x 5x ip'(x) — 3(dx dx)5ip'(x) 

(25) 

= — 5x[5ip(x) — 5xip'(x)] — 3(dx dx)5ip (x). 

Everything in the r.h.s. of equation (25) is well defined at this stage, and if 
we assume that all the terms belong to C3 with /i > 1 [which can be justified 
via Taylor's expansions whenever x € Cf with a > | and <p G C 2 (M)], we can 
conclude that 

3(dxip(x)) = 5xip(x) + 3 (dx dx)ip' (x) 

+ h\3(dxdx)5ip'(x) +5x(5ip(x) - 5xip'(x))], 
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or stated otherwise 

J{dxip(x)) = (Id- Ad)[5xip(x) + J{dxdx)ip'(x)\, 

where we used the fact that 5J(dxdx) = 5x5x to put in evidence the 
fact that we are actually integrating (in the sense of Corollary 2.5) the 
1-increment 5x <p(x) + J{dx dx)(p'(x) which can be thought of as a corrected 
version of the more natural integrand 6x<p(x). It is worth noticing at that 
point that this integral has now to be understood as an integral over the 
(step-2) rough path (x,J~(dxdx)) introduced in [6] and it coincides with the 
notion of integral over a rough path given by Lyons in [17]. 

Remark 2.8. This algorithm has an obvious extension to higher or- 
ders if we assume that a reasonable definition of the iterated integrals 
J{dxdx ■ ■ ■ dx) can be given. To proceed further, however, we need the no- 
tion of geometric rough path (for more details on this notion see [17]) which 
must be exploited crucially to show that some terms are small enough and 
belong to the domain of A. For a more general approach, which does not 
rely on geometric rough-path, see [10]. 

3. Algebraic integration associated to a semigroup. The aim of this sec- 
tion is to set the basis for our future computations: after recalling some basic 
facts about analytic semigroups, we will define a set of increments C* and a 
modified operator 5 adapted to our evolution setting. Then we will give some 
basic calculus rules for (C#, 8) and eventually, we will fix the notation for the 
main application we have chosen, that is the stochastic heat equation. 

3.1. Analytical semigroups. As in [12], we will be able to develop our 
integration theory in the abstract setting of analytical semigroups on Banach 
spaces, whose basic features can be summarized as follows: let (23, | • |) be a 
separable Banach space, and (A, Dom(A)) be a nonbounded linear operator 
on B. We will assume in the sequel that (see [20, Sections 2.5 and 2.6]) A is 
the generator of an analytical semigroup {St',t > 0}, satisfying 

| St | < Me~ xt for some constants M, A > and for all t > 0, 

where | • | also stands for the operator norm on B. Set now A Q = —A. This 
allows us, in particular, to define the fractional powers (A" , Dom(A")) for 
any a£l. 

For a > 0, let B a be the space Dom(A^) with the norm \x\& a = |-A"a;|. 
Since A~ a is continuous, it follows that the norm | • \g a is equivalent to the 
graph norm of A%. If a = 0, then B a = B and A o = Id. If a < 0, let B a be the 
completion of B with respect to \x\s a = \A"x\, which means in particular 
that B a is a larger space than B. We will also set B-oo = U^eR 



16 



M. GUBINELLI AND S. TINDEL 



Among the important facts about these spaces, note the following ones: 
For any a E R and any p > 0, 



(26) Aq P maps B a onto B a+P for all a G R, p > 0, 

(27) |^|Ba — Co,pI x Ib p f° r an x ^ and all a < p. 
Moreover, for all a, j3 £ R, 

(28) y^Af = on £ 7 

with 7 = max{a,/3,a + The semigroup (Sj)t>o also satisfies 

(29) St may be extended to B a for all a < and all t > 0, 

(30) St maps £ Q to B p for all a G R, p > 0, t > 0, 

(31) foralU>0, a>0 \A%S t \ < M a t~ a e- Xt , 

(32) for 0<a;<l, xeB a \S t x - x\ < C a t a \A%x\. 



We will denote with C(B,B') the space of continuous linear operators from 
the Banach space B to the Banach space B' . We let C{B) = C(B,B). In 
order to be coherent with our previous notation, we also set St- S = St s for 
a generic semigroup S, and < s < t < T. 

3.2. Convolutional increments. Let us turn now to the main concern of 
this section, that is the definition of a complex (C*,5) which behaves nicely 
for the definition of our evolution problem. 

Notice that, due to the fact that the operator St^ is well defined only 
for t\ > t2, our integration domains will be of the form S n , where <S n stands 
for the n-simplex 

S n = {(t 1 ,...,t n ):T>t 1 >t 2 >--->t n >0}. 

Let then V be a separable Banach space. The basic family of increments we 
will work with is {C n (V);n > 0}, where C n (V) denotes the space of continu- 
ous functions from S n to V. Observe that an operator 5:C n (V) — >C n+ i(V) 
can be defined just like in (7). In particular, if A E Ci(V) and B £ C2(V), 
the relation (8) is still valid. However, let us see now why 5 is not adapted 
to the resolution of equation (4). 

What made 5 an interesting operator in Section 2 was the simple fact 
that, if F G Cf°(M), then, for t,s(E [0,T] 2 , we have 

(33) [SF] ts = [ f u du with f = F'. 

J s 

However, if St is the semigroup defined at Section 3.1, and if we set 

Ft = f Stufu du for t > 0, / G Cf° (B) , 
Jo 
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then the same kind of relation does not hold true for F. Indeed, for s < t, if 
we define the operator at s ■ B — > B as 

(34) a ts = S ts ~ Id, 

where Id : B — > B is the identity operator, then it is easily seen that 
[5F] ts = F t -F s = a ts F s + / S tu f u du, 

J s 

and hence, in order to get a similar relation to (33) in this new context, one 
should consider an operator 5:C n (B) — >C n+ i(£?), defined by 

[5A] tl ... tn+1 = [5A] tl ... tn+1 - a tl t 2 A t2 ... tn+1 

(35) 

for AeC n (B),{t 1 ---t n+l )eS n+l . 

In the remainder of the paper, we will write 5 A = 5 A — aA, where we made 
use of the convention (16). As in Section 2.1, one can define, for n> 1, 

ZC n {B) = C n (B) n ker(<5) and BC n (B) = C n (B) n Im(5). 

Then the perturbed operator 5 preserves some important properties of 
the original coboundary 5. 

Proposition 3.1. The couple (C*, 5) is an acyclic cochain complex: ZC n +\ = 
BC n for any n > 0. 

Proof. Let us prove first that 5 is a coboundary, that is, 55 = 0. Indeed, 
if F € C n according to the fact that 55 = and thanks to the forthcoming 
Lemma 3.2, we have 

55F = (5 -a) [{5 - a)F] = 55F - 5(aF) - a5F + aaF 

= -5aF + a5F - a5F + aaF = aaF - 5aF. 

Furthermore, it is readily checked that 

(5a)tus = a tu a us , (t,u,s)eS 3 , 

which gives 55 F = 0. 

The fact that Im^,^ = ker^,^ can be proved along the same lines as 

for the (C*, 5) complex [6]: pick A £ C n+ \ such that 5 A = 0, and set B tl ...t n = 
^i...t n s, with s = 0. Then 

[5B] tl ... tn+1 = [5A) tl ... tn+ls + (-l) n+1 A tl ... tn+1 - a tl t 2 At 2 ...t n s 

= [5A] tl ... tn+lS + A tl ... tn+1 = (-ir+ 1 A h ... tn+1 . 
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Thus, setting C = (-l) n+1 B, we get 5C = A. □ 

The cochain complex (C*,<5) will be the structure at the base of all the 
constructions in this paper. Let us also mention at this point that, when 
the meaning is obvious, we will transpose the notation of Section 2 to our 
infinite-dimensional setting. Furthermore, whenever this does not lead to an 
ambiguous situation, we will write C n instead of C n (B). 

Let us give now a simple and useful extension of Proposition 2.6, which 
has already been used in the last proposition. 

Lemma 3.2. Let L G C n _i(23) and M G C 2 (£(£)) . Then 

5(ML) = 5ML - M5L. 

Proof. Let G tl ... tn = M tlt2 L t2 ... tn . Then 

n+1 

[SG] tl ... tn+1 = X)( _1 ) lG! tvti-Wi 
i=l 

n+1 

= -M t2 t 3 L t3 ...t n+1 + M tlt3 L t3 ... tn+1 + J2(-lfM tlt2 L t2 ^ tn+i 

i=3 

n+1 

= [SM] tlt2t3 L t3 ... tn+1 + M tlt2 Ysi-^K-ii-tn+v 
which yields our claim. □ 

3.3. Computations in C*. Here again, like in Section 2, we will try to 
move from a smooth setting to an irregular one. And we will start by giving 
the equivalent, in our new setting, of Proposition 2.3, which will require first 
the introduction of some analytical structures on the spaces C n . 

3.3.1. Holder type spaces. First of all, we have to define some Holder 
type subspaces of Ck, k < 3, related to the spaces B a , a € R: for fj, > and 
g G C2(B a ), we set 

(36) l|ff|L,a = sup jf ts g " and C% ,a = {g G C 2 (B a ); \\g\\u,, a < oo}, 

and the definition above also induces some seminorms on C\\ for 7 > 0, a G M, 
we say that / G Cj' a if 

II/II 7 ,q = II^/IIt.o < °°- 

Another useful subspace of C\ will be C°' Q , the space of bounded paths in 
B a with the supremum norm \\f\\o, a = sup te[0)T ] \ft\B a - 
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As far as C3 is concerned, C^ a can be denned in the following way: set 

\htus\j3 a 



t,u,s€S 3 \t-u\^\u-s\P 

(37) 

i*,<* = inf I X] ll^llw,M-w,a! h = Y^hi,0<pi<n\, 

i i 

where the last infimum is taken over all sequences {hi}i such that h = ^ hi 
and for all choices of the numbers pi <G (0, p). Then || • || M;Q is again easily 
seen to be a norm, and we set 

C%' a = {heC 3 (B a );\\h\\ lx , a <oc}. 

In order to avoid ambiguities, we shall also denote in the sequel by J\T[f;Cj] 

the K-Holder norm on the space Cj, for j = 1, 2, 3. For £ € Ci(V), we also set 

AA[C;C 1 (y)]=su Pse[0;T] ||C s ||y. 

Eventually, we will need to introduce a slight extension of the spaces we 
have just defined above: for j = 1,2, let £j' a be defined by 

(38) £f Q = p| q~ £ ' a+£ , 

£</iAl~ 

where £</iAT stands for the condition e £ [0 ;/ u] n [0, 1), and where the 
intersection is considered along any arbitrary family {0 < e± < ■ ■ ■ < e n < p A 
l - } for n > 1. Obviously, some families of operators will play an important 
role in the sequel, and this will lead us to the following specific definitions 
for operator-valued increments. 

Definition 3.3. For p > and a,/3 G R, we will call C^C?^ the space 
C%{£(Bp;B a )), and will denote by the space 

£</lAI~ 

where the intersection is still considered along any arbitrary finite family 
{0 < E\ < ■ ■ ■ < e n < p A l - } for n > 1. The natural norm on &CP^ will be 
defined by 

(39) \\A\U, a = sup M% 

and when we consider some Hilbert-Schmidt operators, the corresponding 
spaces will be denoted by C%£^ and S^C^ . 
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3.3.2. The convolution sewing map and related properties. Here is a first 
proposition showing how the analytical structures introduced above interact 
with our previous algebraic notation. 

PROPOSITION 3.4. If fi > 1, then for any a£l, ZC%' a = {0}. 

Proof. Take h £ ZC^ 01 . Then, according to Proposition 3.1, there exists 
/ £ C\ such that h = 5f. Consider the telescopic sum 

n 

h ts = (Sf) ts = S tti+1 (Sf) ti+lti , 

i=0 

with respect to the partition Hf s = {to<«<n+i :to = s, t n+ \ = t} of the interval 
[s,t]. Since 5f £ ZCi^' 01 with fi > 1, we have 

n n 

\(sf) ts \ Ba <^2\(sf) ti+lti \ Ba < \\sfWw 

i=0 i=0 

which converges to zero as the size of the partition goes to zero. Since t, s 
are arbitrary, we have 5f = h = in Cr, ,a . □ 

We can now state and prove the equivalent of Proposition 2.3 in our 
evolution setting, which is the main aim of this section. 

Theorem 3.5. Let fj, > 1, a € M. There exists a unique sewing map 
A : ZC^' a — > ° su °h that 5 A = Id^ 3 . Furthermore, for any < e < /i A 1~, 
there exists a strictly positive constant c^ :£ such that 

(40) \\kh\\^ £)a+£ <c^ £ \\h\\^ 

for any h £ ZC^' a . 

Proof. Like in the proof of Proposition 2.3, we will divide our compu- 
tations in two steps below. 

Step 1 : The uniqueness part of our theorem simply stems from the fact 
that if we have 5a = h and 5a' = h with a, a' £ C 1 ^ ' a , then b = a — a' £ ZC^ a 
and since fj, > 1, by Proposition 3.4, we must have b = 0. 

Step 2: The existence part can be adapted from Proposition 2.3, and 
we will construct a process M £ ,a such that 5M = h starting from any 
B € C2{B a ) satisfying 5B = h (this increment B exists thanks to Lemma 
3.1). Now, similar to (13), we will set, for a given n > 1, and (t,s) £ £2, 

2 n -l 

M U =B ts - S tr ? +1 B r n +i:r n, 
i = 
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where s,t and r" have been defined at (12). Then Af& = and 
Af2 +1 -AfS 

2 n -l 

= / (^r- n+1 B n+1 n+1 — S, n+1 B n+1 n+1 — S+ n+1 B n+1 n+1 ) 
Z — / tr 2i+2 '2i+2> r 2i "2i+2 r 2i+2>'2i+l w 2i+l r 2i+l>'2s 

i=0 
2 n -l 

= / S, n + 1 (B n + 1 n+1 — B n + 1 n + 1 ~ B n + 1 n+1 ) 

Z— ' L, 2i+2 '2i+2''2i r 2i+2>'2i+l r 2i+l>'2i 
1=0 

— 5, n+1 [S n+1 n+1 — Id] B n+1 n+1 . 

11 1i+2 '2i+2'2i + l ' 2i + l'' 2i 

Thus, according to the definition (35) of 5, we get 

2 n -l 

M^ +1 — M™ = 5", n + 1 [(55) n+1 n+1 n + 1 — a n+1 n+1 -B n + 1 n+l] 

Z — ' Ir 2i+2 r 2i+2' T 2i+l> r 2i r 2i+2> T 2i+l r 2i+l' T 2i 

i=0 



2™-l 2™-l 



n + 1 „n+l n+1 
)' 2i_l_1 i' o 



~~ 7 , + 1 {SB) n + 1 n+1 n+1 — > S, n + 1 /l 

* — ' u 2i+2 f 2i+2>'2i+l>'2i — ' t/ 2i+2 

i=0 i=0 

Hence, for any e < //, we get, invoking (31), 

2 n -l 

K +£ (M« +1 -M«)|<c e 53 It-r^riAl^lt-^r 

8=0 

r 1/ — sl^~ e l/)l f 1 
2»C**-i) 7o 

which gives, like in Proposition 2.3, that Mt s = hnin^oo M^ s exists, and is 
an element of Now, the fact that 5M = h can be shown analogously to 
the case of Proposition 2.3, and the proof of (40) is straightforward. □ 

A direct consequence of the existence of the A-map is a result of conver- 
gence of finite sums. 

Corollary 3.6. Let g € C2 such that 5g £ C^' a for some fj, > 1. Then 
the 1-increment 5f = (Id — A5)g € satisfies 

n 

(Sf)ts = lim y^Stti+xgti+xti, 
|iits| ^° i=0 

for all (t, s) G ^2- 



Proof. It follows the lines of the proof of Corollary 2.5. □ 
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We will now define an equivalent of the iterated integrals of Section 2.2 
in our convolution context: consider some smooth functions g £ C^°(C(B a )) 
and / € Cf°(B ), for some q£K. Then J(dg f) will be defined as an element 
of 62° (B a ) by Jts{dg f) = f^dg v f v , for (t,s) € 52. We will also need some 
integrals of processes weighted by the semigroup S, defined as follows, for 
< s < t < T: 

Jts (dg f)= s tv dg v f v . 

J s 

Once these elementary blocks have been defined, the iterated integrals 

(41) J{d*-f n ---d^f x ) for /„,..., f 2 eC?(£(B a )),fieC?(B a ), 

where d* j fj stands for any of the increments of the form dfj or dfj, can be 
defined recursively along the same lines as in Section 2.2. In particular, the 
operator-valued increment J(dgS) is defined by 

t 



Jt s (dgS) = / S tu dg u S; 

J s 



The relations between 5 and these integrals, which will be useful for our 
purposes, can be summarized in the following: 

Proposition 3.7. Let aeR, and geC^°(C(B a )), f eC^°(B a ). Then 
5(J(df)) = 0, 5(J(dgf)) = 0, 5(J(dg5f)) = J(dgS)5(f) 

and 

5(J(dgdf)) = J(dgS)J(df), 5(J(dgdf))=J(dg)J(df). 

Proof. The proof of these results is elementary. We will give some 
details about the last relation for sake of completeness. For any (t, u, s) 
invoking the definition of 5, we have 

[5(J(dgdf))] tus = [5(J(dgdf ))] tus - a tu J us (dgdf) 



t rv pt rv 

Stv dg v / df w - S tv dg v \ df w 

J s J u J u 

/u pv pu pv 

S U v dg v / df w - [S tu - H] / S uv dg v / df u 
J s Js J s 



= Jtu(dg)J us (df), 
which proves the claim. □ 
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3.4. Fractional heat equation setting. In this section, we will give the 
general setting under which we will try to define and solve the stochastic 
heat equation driven by an infinite-dimensional fractional Brownian motion: 
as mentioned in the Introduction, the main application we have in mind is 
the situation where A = A — Id, and A is the Laplace operator on the circle S, 
assimilated to [0, 1]. This operator can be diagonalized in the trigonometric 
basis of L 2 ([0, 1];C), namely {e n ;n G Z}, where 

e n {x) = e 2mnx , xe[0,l]- 

Associated to these eigenfunctions are the eigenvalues \ n = — 1 — (27rn) 2 . We 
have chosen to deal with A = A — Id instead of A itself for computational 
convenience, since this choice avoids the problem of a null eigenvalue for con- 
stant functions. Notice that in this case, A is the generator of an analytical 
semigroup, and all the constructions of Section 3.1 goes through. Then B a 
can be identified with H a , the usual Sobolev space based on L 2 ([0, 1]), for 
the definition of which we refer to Adams [1], and {St',t > 0} stands for the 
heat semigroup, which admits a kernel Gt{£,,rj) for t > and £,77 G [0, 1]. In 
this context, set Gf (£,77) for the kernel of the operator A®St, and G"(£,rj) 
for the kernel of the operator Aq . Then, for a£l and j3 > 0, Gf and G@ 
admit the following spectral decomposition: 

Gf(C,7 1 ) = Y,Ke- tXn e n (C)e n (r ] ) and 

(42) 

G /3 (e,r ? ) = ^A-^e n (e)en(r/). 

Let us specify now the noise X we will consider: we will try to stick to the 
existing literature on the topic, and choose a fractional Brownian noise in 
time, defined on a certain complete probability space (£1, J-, P), which will 
be homogeneous in space, with a spatial covariance function Q. Namely, X 
will be a centered Gaussian field indexed by functions on [0,T] x [0, 1], such 
that if 4> and ip are smooth enough, then 

E[x(<f>)xy,)] 

(43) . , f s 

= ch I / Q(£ - '?)0(«,C)V ; ( v > r /) d^drj I \u - v\ 2H 2 dudv, 

J[0,T] 2 \J[0,1} 2 J 

with ch = H(2H — 1), for H > ^. Notice that, in order to simplify our 
statements, we will generally assume that Q can be decomposed itself on 
the basis {e n ;n G Z} in the following way: 

(44) Q(0 = ^ q n e n {i) with q n = X~ u , for v G [0, 1), 
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and notice that the case v = corresponds to the white noise in space, while 
the case v > 1/2 corresponds to a noise admitting a density in space. Some 
explicit construction of such kind of noise, as well as an account on the 
related stochastic calculus, can be found in [24]. The methodology we will 
develop in the rough case will also enable us to handle the Brownian motion 
case, which means a covariance structure given by 



(45) E[X(4>)X{il>)]= [ <t>(u,OQ(Z-riMu,n)dtdn\d u . 

J[0,T] \J[0,1] 2 J 

We give here a slight extension of a result result of [6], which will be used 
below to prove existence of regular versions of some stochastic processes, 
following the well-known approach of Garsia-Rodemich-Rumsey. The proof 
is conceptually similar to that appearing in [6] but there is a small technical 
difficulty due to the fact that convolutional increments are one-sided and 
which forces us to follow the scheme of the proof of the GRR inequality in 
Stroock's book rather that which can be found in [6]. 

In order to state this extension, we shall introduce for the first time a 
variant of the operator 5, called <5, acting on operator-valued increments 
which turns out to be useful in the sequel, and which is defined by 

(46) 5Q = 5Q-Qa = 5Q-aQ-Qa, QeC*(C(B)). 

With this additional notation in hand, our regularity lemma is the following 
below. 

Lemma 3.8. For any 7 > 0, a,/3 E 1R and p>l, there exists a constant 
C such that for any R € C.2{C^ ,a ), we have 

( 47 ) ||-R|l7,/3,a < C(^ 7 +2/p,p,/3,a(- R ) + P-Rl^.a), 

where 



s 2 



iRtsi ^ a ydtd S 



t-sp 



l/p 



Proof. As in [6], this result is a direct consequence of a more general 
Lemma 3.9 below by choosing <&(x) = x p and p(t) = t 7+2 / p . □ 

Lemma 3.9. Letp and ^ be strictly increasing, continuous functions on 
R + satisfying p(0) = ^(0) = and ^/(x) — > 00 as x — > 00. Then there exist a 
constant K such that for any a,/?El and any ReC 2 {C^ a ) for which 

u= ff J\^\m ) < n (h< x 

J J0<s<t<T \p(t-S) 
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and 

sup ISRtusl^K^f 4C )p{t-s), 0<s<t<T 

s<u<t \{t-sYJ 

for some constant C < oo, then 

\R ts \i3,a<8K 1^ \~i(^!pjp(du) + 9K J* * ^(^pjp(du) 

for allO<s<t< T. 

Proof. The proof follows closely Stroock's proof of the Garsia-Rodemich- 
Rumsey inequality. First, show the estimate for T = t = 1 and s = 0. For a 
sequence of times t n ,s n such that t > t^+i > ifc, s < Sk+i < Sk and to = so, 
we have 

Rts = Rtt St s + Su Rt s + SRtt Q s, 

n n 

Rtt = R tt n+1 St n+1 t + ^ s tt k+1 Rt k+1 t k St k t + ^ S R tt k+1 t k St k t 

k=0 k=0 

and 

n n 

Rt s = St 0Sn+1 R Sn+ls + Stos k Rs k s k+1 Ss k+1 s + ^t s k ^Rs k s k+1 s- 

k=0 k=0 

Next, we choose these times as follows. Let I(v) = J Q V ip( ^yz^y ) du. For 

any s n , define d n by the equation 2p{d n ) = p(s n ). Remark that since I(t) dt = 
U there exists to sucn that I(to) < U . We claim that there exists s n+ i € 
(0, d n ) such that both inequalities 



dn \P(Sn ^n+l) / d n 

hold. This is always possible since, if we call A n C (0, d n ) (resp. B n ) the set 
of s n+ i where the first (resp. the second) fail, we have 

f 2U 

U> ds n+1 I(s n+1 ) > — \A n \ and 

I{s n )>( ds n+1 *( \f^\ )> 2 ^\B n \ 

so we must have \ A n \ < d n /2 and \B n \ < d n /2 which means that (0, d n )\(A n U 
B n ) has positive measure. Then since 

p(s n - s n+ i) <p(s n ) = 2p(d n ) = 4(p(d n ) -p{d n )/2) < 4(p(d n ) - p(d n+1 )) , 
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we have 



and 



\R Sn s n+1 \p,a < ^ {^^-)p( S n ~ *n+l) 



then we have 

oo oo 
\R t0 s\p,a<M 2 J2\R \ + Mj2\SR, 



n=0 n=0 



n=0 Jd n+1 V « / n=0 V " / 

where we used the fact that there exists M > 1 such that |5 t |q, jQ < M for 
any a and any r > 0. 
Similarly, we find 

\Rtt \p,a<4M 2 ^{ — \p(du)+AMj^ <S>-U — \p(du). 

So using that 

(4C \ f t ~ s /4C\ 

we obtain 

\Rts\p,a<8M 2 _jp(du) + 9M jf 

Now it is not difficult to extend this to generic < s < t < T. □ 
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4. Young theory. We are now ready to analyze the Young integration in 
the evolution setting along the same lines as in Section 2.3: we will first define 
the integral J(dxz) for two Young paths x,z in an abstract setting. Then 
we will solve Young SPDEs, and eventually, check our main assumptions in 
the fractional heat equation setting of Section 3.4. 

4.1. Young integration. The extension of the notion of integral weighted 
by an analytical semigroup will be performed through the following algo- 
rithm, which will be used in fact throughout the remainder of the paper: 

(1) Assume first that x is a regular operator- valued increment, and z a 
regular B- valued function and let Jt s {dx z) = J St u dx u z u , for (t,s) € 
<S2, as an element of 

(2) Through the application of 5 and A, try to get an expression for J{dxz) 
which depends only on minimal regularity requirements for x and z. 

(3) Extend the notion of integral using the previous step, and see that it 
induces the convergence of some well-chosen Riemann sums. 

Here is how this general strategy can be implemented here: suppose for 
the moment that x is a smooth operator-valued function and z a smooth 
function. Then it is easily checked that 

(48) J{dx z) = J{dx S)z + J{dx 5z). 

Note that in this last equality appears for the first time an incremental op- 
erator which will play a fundamental role in the sequel, namely the operator 
X 1 eC 2 (C(B)) defined by 

(49) Xl = J ts (dxS) = j S tu dx u S us . 

And here is an important point of our strategy: the noise x does not appear 
by itself but always inside a convolution of the form (49), so its action is 
milded by the regularizing properties of the semigroup. 

Applying 5 to the last term of equation (48) and invoking Proposition 
3.7, we get 

5[J{dx5z)]=J{dxS) 5z = X 1 5z. 

If the 2-increment X 1 5z is small enough, namely if X 1 5z E ' for some 9 
and some \i > 1, then we can express J(dx z) as 

(50) J(dx z) = X 1 z + AfX 1 6z] = (Id - M)[X 1 z\. 

The last equality is justified by noting that when x is a smooth incremental 
operator, we have 5X 1 = X 1 a (i.e. 5X 1 = 0), and thus by Lemma 3.2 one 
obtains that 5(X 1 z) = —X 1 5z. 
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Let us turn now to the second point of our general strategy, which consists 
in inverting the process which leads to (50): indeed, if we can define properly 
the right-hand side of (50), then we will be able to extend the notion of 
integral by a procedure which is coherent with the basic properties required 
to any integral J . Notice that this step only relies on the definition of an 
operator X 1 associated to x, satisfying 5X 1 = 0, and such that X 1 is regular 
enough. This will be formalized in the following theorem (recall that the 
space B-oo has been defined at Section 3.1). 

Theorem 4.1. Let then x be a path from [0,T] to B-oo such that the 
operator X 1 associated to x is well defined as an element ofC^C 13 ' , where 
/3, k are positive constants, and a G R. We also assume that X 1 satisfies the 
algebraic relation 5X 1 = X 1 a. Let z G &y , with k + r\ > 1, and set 

(51) J{dxz)=X 1 z + k[X 1 5z] = (Id - AS^z}. 
Then 

(1) J(dxz) is well defined as an element of ■ 

(2) For a constant c > 0, we have 

\\J{dxz)\\^ a < HX^I^adlzHo^ + ^1^11^), 

where the norm \\ ■ || Kj /3 jQ . has been defined at relation (39). 

(3) It holds that, for any < s < t <T 

n 

JUdxz) = Jim VSh+i^+iA^' 

where the limit is over all partitions Ut s = {to = t, . . . ,t n = s} of [s, t] as 
the mesh of the partition goes to zero. 

Proof. Since X l z is a well defined element of C^' 01 , in order to show 
that the r.h.s. of equation (51) is well defined, it only remains to check 
that X 1 5z is in the domain of A. However, since we have assumed that 
X 1 G C^C 13 ' and 5z e , we obviously get that X l 5z G C^ +V ' a . Thus, 
according to Theorem 3.5, X 1 5z G Dom(A), yielding the first assertion of 
our theorem. 

Moreover, thanks to the second part of Theorem 3.5, we have 

(52) IIApT^H^cjA: 1 !! 
and it is also readily checked that 

(53) HA^zll^a < HA" 1 || K ,/3, Q ||£||o,/3, 

which shows our second claim, by using equation (52) and equation (53) to 
estimate the r.h.s. of (51). Eventually, the third part of the theorem is a 
direct consequence of Corollary 3.6. □ 
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Remark 4.2. It is worth stressing at this point some elementary prop- 
erties enjoyed by the extension of the notion of integral given by Theorem 
4.1: 

• The third part of the theorem states that Jt s [dxz) is associated to some 
natural Riemann sums involving the processes x (through X 1 ) and z. 

• The arguments leading to relation (50) also show that, in case of some 
smooth processes x and z, our integral J ts (dxz) coincides with the usual 
one. 

These first properties seem to imply that our integral extension is a reason- 
able one. 

4.2. Young SPDEs. Recall that we wish to solve an equation of the form 

(54) dy t = Ay t dt + dx t f(y t ), te[0,T], 

with an initial condition yo = ip £ B K , where x is an operator- valued process 
which represents our noise and / :B — > B is a (possibly) nonlinear regular 
map. As mentioned in the Introduction, we will consider equation (54) in 
the mild sense, that is, we will say that y is a solution to (54) if, for a given 
k > (specified below) we have y € C*' K and if, for any t S [0,T], yt satisfies 

(55) yt = S t ip+ S tu dx u f(y u ) = Styo + Jto(dxf(y)), 

Jo 

where the integral Jto{dx f{y)) is understood in the sense of Theorem 4.1. 
In fact, we will focus here on a slight extension of the problem given by (55): 
we will search for a (unique) process y £ C^' K satisfying, for any (t,s) € S2, 

(56) yt = S ts y s + Jts(dxf(y)), Vo = i>, 

from which one recovers obviously (55) by taking s = 0. Now, (56) can be 
expressed in terms of convolution increments, since it is equivalent to the 
following one: 

[5y}ts = Jts(dxf(y)) = [(U-A6)[X 1 f(y)]} ts for {t,s)eS 2 and 

(57) 

yo = i/>, 

which sticks better to the algebraic formalism introduced in the previous 
sections. 

Let us specify also some of the assumptions under which our computations 
will be performed: first of all, the incremental operator X 1 defined by (49) 
will be assumed to be in the following class. 
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Hypothesis 1. Assume that X 1 £ nC£°£ K ' K for some 7 > k > 

k > 1/4 such that 

7 + k>1, 7 — k>kq, k<1/2. 

Notice that in the hypothesis above, the condition k < 1/2 is somehow 
redundant. Indeed, if 7 > k + kq > 2k, this forces the relation k < 1/2. 

As far as the function / is concerned, we will also assume that the fol- 
lowing holds true. 

Hypothesis 2. Let k be the strictly positive constant defined at Hy- 
pothesis 1. We assume that the function f : B K — > B K is locally Lipschitz, and 
satisfies \f{x)\s K < c/(l + |£|e K ). Furthermore, we suppose that f can also 
be seen as a map from B to B, and when considered as such, it holds that f 
is globally Lipschitz. 

With these assumptions and notation in mind, we are now able to solve 
our evolution equation in the Young sense. 

Theorem 4.3. Assume Hypotheses 1 and 2 hold true, and that ip € B K . 
Let Cl' K be the subspace of C\ defined by the norm 

(58) || ^11*, k — ||z||o,k "t~ I \°~ z || K, re- 

Then there exists a unique global solution to (57) in C^ K . Furthermore, this 
solution enjoys the following properties: 

(a) For any t £ [0,T], y t can be written as yt = St*p + (Sy)to- 

(b) Let us call $ the map (ip^X 1 ) 1— > ^(^X 1 ) = y, where y is the solution 
to (57). Then $ is Lipschitz continuous from B Ko x (C^C°~ K nC^C*'* 1 ) 
to C 1 . 

Proof. A classical fixed point argument will be sufficient to obtain 
the global solution. Let us introduce the map T:C\' K — > C\ ,K defined in the 
following way: if y € we set T(y) = z, where z satisfies 

[Sz}ts = Jts(dxf(y)) 

(59) =X 1 f(y) + [A[X 1 of(y)]} ts for (t,s)eS 2 and 
z =ip. 

Let also B be the ball defined by 

(60) B = {y;y = if,, ||y||, >(t < 2(1 + |^| B J}. 

Then the fixed point argument can be decomposed into two usual steps: 
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(1) Show that, on a small enough interval [0,T], the ball B is left invariant 
byr. 

(2) Prove that T, restricted to the ball B, is a contraction. 

We will mainly focus, in this proof, on the first of these steps, since it contains 
most of the technical difficulties associated to our method. 

Take y € B, and let us show that z = T(y) € B whenever T is small enough 
(recall that S2 depends on the parameter T). To this purpose, we will first 
bound the term A[X 1 <5/(y)] in (59). Recall that 

(61) [5f(y)] ts = [5f(y)] ts -a ts f(y s ), 

and let us estimate the terms in the right-hand side of (61) separately: on 
one hand, recalling the notation of Section 3.3.1, and thanks to the fact that 
/ is Lipschitz on B, we have 

\[8f(y)]ts\B < cf\[6y]ts\B < cf(\[5y]ts\B + \(HsVb\b) 

(62) <c f [\\y\\ K , + \\y\\ 0iK }\t-s\ K 

< Cf\\y\\* tK \t - s\ K , 

where ct is a positive constant which may change from line to line, but 
which depends only on /. On the other hand, according to Hypothesis 2, it 
is readily checked that \f(y s )\l3 K < Cf(l + ||y||o,/t)- Thus, invoking (32), we 
obtain that 

(63) |ot./(y.)lB<c/(l + ||l/lk«)|t-sr, 
and plugging (63) and (62) into (61), we get 

(64) \[Sf(y)}t s \B<cf(l+\\y\U, K )\t-s\ K . 

However, we know that X 1 € and this fact, together with the last 

estimate, yields 

\[X x Sf{y)\ta.\B. K < c f \\X% fi! _ K {l + \\y\U, K )\t - up\u - s\ K . 

Furthermore, by Hypothesis 1, we have 7 + n > 1. This means that Theorem 
3.5 can be applied here to obtain that A(X 1 5f(y)) € £2 +K ~ K - in particular, 
invoking the definition (38) of the space £2 +K '~~ K , and since 2k < 1 and 
k < 7 - k, we get MX 1 6f) G C%~ K ' K C C^ K . Moreover, 

(65) ||A(A' 1 ^/)|| lto ,«< C/ ^ /s ||A' 1 ||^o,_«(l + ||y||*,«). 

A bound similar to equation (65) can be found for the term X 1 f(y) ap- 
pearing in the definition of 6z in equation (59). Indeed, owing to the fact 
that X 1 £ C2°C K,K and that / has linear growth in B K , we get 

\xlf(y s )\B K < c f \\X% 0iKtK (l + \y s \ B J\t - s p 

(66) 

<c / ||X 1 |U 0A)K (l + ||y||* iK )|t-s| K0 . 
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Hence, plugging (66) and (65) into (59), one obtains that ||<5z|| KOi/t < ct x 1 7 «(1 + 
||y||* «)• Note here a crucial point: starting from y € C^' K , we have constructed 
z € C*°' K with e = kq — k> 0. This little regularity gain can be used in order 
to write 



(67) \\6z\\ K)K <c ftX x >%K (l + \\y\\* )K )T s . 



Now, the quantity T £ can be made arbitrarily small as T — > 0. Moreover, 
recall that we still need a bound on H^ll*^ defined by (58), and thus an 
estimate on ||z||ok is needed at this point. However, it is easily checked that 



which yields that, whenever cT £ < 1/2, the ball B defined by (60) is left 
invariant by the map V. 

Now that the invariance of B has been shown, the contraction property 
for r in a small interval [0,T] is a matter of standard arguments, and is 
left to the reader for sake of conciseness. Let us just mention that / is only 
supposed to be locally Lipschitz when considered as a function from B K to 
B K . However, we are able to establish the contraction property here, due to 
the fact that we are confined to the ball B. This gives the existence and 
uniqueness result for equation (57) in the small interval [0,T] whose size 
does not depend on the initial condition ip. The construction of a global 
unique solution from the solution in [0,T] is also quite standard, and its 
proof will be omitted here. □ 

4.3. Application: the fractional heat equation. Let us see now how the 
abstract results of Section 4.2 can be applied in the case of the heat equation 
driven by a fractional Brownian motion defined at Section 3.4. Recall that 
this means that we wish to solve equation (55) in case A = A — Id, where 
A is the Laplace operator on the circle, x is a fractional Brownian motion 
defined by the covariance function (43), B K stands for the usual Sobolev 



space on [0, 1], and / ■ B K y 13 K is defined by [f(y)](0 = a(y(0) for £ € [0, 1] 



and a smooth function a :M— >R. In other words, we will try to solve the 



(69) y(t,0= G t (C,v)^(v)dri+ / G t - s (Cv)X(ds,dn)a(y s ( V )), 



where the last integral has to be understood in the sense of Theorem 4.1. 
Notice that we have chosen here a multiparametric formulation for our equa- 
tion, for computational purposes. However, as mentioned in the Introduc- 
tion, this setting can be translated easily into the infinite-dimensional one. 



(68) \zt\ BK < \SMb k + \(Sz) t0 \B K < Ms K +T Ko \\5z\\ KOiK 

Putting together (67) and (68), we finally get, on [0,T], that 
\\z\\*,k < \iP\b k + c(l + \\y\\*,K,)T £ with c = c f}X ifi 



equation 
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Now, the application of Theorem 4.3 in this context amounts to define an 
incremental operator X 1 related to our problem, and then to show that 
Hypotheses 1 and 2 are fulfilled. 

Let us give then a natural definition of the operator X 1 associated to our 
equation: we will set, for ip G B and (t, s) G 1S2, 

[xim) = [Jts(dxs)]m 

(70) 

which has to be understood now in the Wiener sense, as a centered Gaus- 
sian random variable whose variance is given by (43). In this context, the 
regularity result we obtain on X 1 is the following. 

Proposition 4.4. Let X be an infinite- dimensional fractional Brown- 
ian motion defined by the covariance function (43) for a given H > 1/2, with 
Q given by (44) for v G [0, 1). Suppose that H + v/2 > 3/4, with the conven- 
tion v = v A (1/2). Let X 1 be the incremental operator defined by (70). Then 
for any 7 < H — 1/4 + v/2, k G (1/4, 1/2), kq = 7 — k and 7 < H we have 



almost surely. 



Remark 4.5. The reader will probably notice that the assumption k > 
1/4 is not necessary in order to prove the proposition above. However, we 
include it already at this stage, since this restriction is crucial for Proposition 
4.10 to hold true. 

The proof of Proposition 4.4 relies on the following elementary lemmas, 
that we label for further use. 

Lemma 4.6. For any a < (3, such that a + f3 > 1/2, there exists a con- 
stant C such that 

i,j:i+j=k 

where f3 = min(/3; 1 /2) . 

Lemma 4.7. Let a and b be two positive constants, and H > 1/2. Then 
the integral 

-1 rl 

/ \u-v\ 2H ~ 2 \2-u-v\~ a \u + v\~ b dudv 
'0 Jo 

is finite whenever 2H — a > and 2H — b>0. 
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We leave the easy proof of these results to the reader. 

Proof of Proposition 4.4. We need to prove that the r.v. X 1 has a 
version with the claimed regularity. For random operators, up to our knowl- 
edge, no standard method is available to prove regularity properties. So we 
have chosen the following simple (though arguably nonoptimal) strategy 
in order to obtain a regular version: first, we determine the kernel asso- 
ciated to the operator X , then using the kernel we estimate its Hilbert- 
Schmidt norm in some L 2 space. This will be enough to apply the modified 
Garsia-Rodemich-Rumsey Lemma 3.8 and conclude the proof. We will de- 
velop now this strategy into several steps, discussing in detail the proof of 
X 1 € Cj/3yg K . The other pathwise statements can be proven similarly. 

Step 1: Definition of a random kernel. For (t,s) G £2, is consid- 
ered as an operator from B = L 2 ([0, 1]) to B~ K , and thus \\Xf S \\hs,b^-B- k = 
\\A~ K X^IIhs^-s-B) which is the expression we are going to evaluate. Pick 
ip € B smooth enough. Applying Fubini's theorem for the fractional Brown- 
ian motion, we get 



J 

where the kernel GJ* U has been defined at Section 3.4, and where K ts {^,rj) 
is the random kernel on [0, l] 2 defined by the Wiener integral 



Js Jo 

Hence, the Hubert-Schmidt norm of X^ s , seen as an operator from B to 
£>_ K , will be given by 




x 







(71) 
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= c H [ [ ( [ G^z)G u . s ( V ,z)Q(z-z) 

Js Js \J[0,1] 2 

x Gt* v (£,z)G v - s (r],z)dzdZ\ \u - v\ 2H ~ 2 dudv. 

Furthermore, for z,z € [0, 1], it holds that 
l 

Gt-u\£> z )Gt-v(£i *0 d£ = G 2t * u _ v {z, z), 

o 
l 

G u - s (7],z)G v - s (r],z)dr] = G u+v - 2s (z, z). 

o 

Thus, going back to relation (71), we obtain 
A ts = ^[ll^sllls] 

= c H / ( / Q{z - z)G u+v -2s{z,z)G2t- u - v {z,z)dzdz 

Js Js \J[Q,1] 2 

x |n - v\ 2H ~ 2 dudv 

\2H-2 



= ch / F(u,v)\u — v\ dudv, 
Jo Jo 

where we have set e = t — s, and with F : [0, e} 2 — > M+ defined by 

(73) F{u,v)= [ Q(z-z)G u+v (z,z)G^ u _ v (z,z)dzdz. 

J[o,i] 2 

Furthermore, plugging the definitions (42) and (44) into (73), and invoking 
the fact that {e n ;n £ Z} is an orthonormal basis of L 2 ([0, 1]), we get 

F(u,v)= Yj \- v \- 2K e- x ^ u+v h- Xli2£ - u - v \ 

m,n,l£D 

where D = {m, n,l6Z 3 :m + n + I = 0}. Then 



A ts = C H }2 X n U K 2K / | 7t _„|2-2ff dudV - 



m,n,l£D 



Owing now to the fact that x i— > x a e x is a bounded function on IR + for any 
a > 0, we obtain, for a constant c which may change from line to line, 

A ts < c £ \-»k 2k k:[ f dudv 

JO JO 

(74) 

< 2H-a \^ \—u\-2K\—a 

— c£ z_> n i m ' 

m+n+l=0 



Jo \u — v\ 2 2H (u + v) a 
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where we have used Lemma 4.7 under the condition a < 2H. Let us now 
analyze the sum. Of course, we can write 

E\ — v \ — 2k \ — a ^ \ r \— 2k \ , \— a\— v 
A n A l A m S I 2_> m n • 

m+n+l=0 l,k:l+k=0 m,n:m+n=k 

Moreover, taking a = 1/2 — v + r) for some small rj > and using Lemma 4.6, 
we have 

E W"v<c E A r -Ar p+1/2 = c e W. 

m+n+Z=0 Z,fc:Z+fc=0 Z,fc:Z+fc=0 

and this sum is always finite under the condition k > 1/4. Then, going back 
to (74), we have found that A ts < ce 2 ^' , for any f = H-a/2 < H- 1/4+9/2, 
where we recall that v = inf(z/; 1/2). 

5*iep 3: LP estimates. We will prove now that, for any p > 1, we have 

(75) £[11^11^] < - s) 2 ^ for < S < t < T. 

Indeed, a simple application of Holder's inequality yields 



E[\\Xl\\*]= I E 



[0,1] 



.j=i 



< / ]lE 1 '*[K%(€ i ,T H )]dt 1 d m ---d$ p dri P , 

J\0,1] 2 p ~T 



'[o,i] 2 - i= i 

and since Kt s (rji,r]i) is a Gaussian variable, we get 



•Ao,i] 2 "f = i 



E[KKS, V )]dSd V ' 
[o,i] y 

= <^[p&||U 
which easily yields (75). 

Step 4 '■ Conclusion. Recall that X 1 is considered as an element of ^(-Cgg K ). 
We can use now inequality (47), which can be read here as 

(76) \\X l \\ 

7,0,— k +2/p,p,0,- 

.(X 1 ) + \\SX^\ 

7,0,— k] j 

in order to bound WX 1 || 7 ,o,-k for any j <j' < H — 1/4 + u/2. Indeed, if p is 
large enough, we have that 7 + 2/p < 7', and the term U^+2/ P)P q K (X ) is 
easily handled thanks to (75). This yields 

(77) Ep^^^iX 1 )] < 00. 
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We are now left with the estimation of ||<5X ||^. However, remember that 
5X 1 = in case of a regular signal x, and it is readily checked that this 
relation is still valid in the current fractional Brownian setting, so this term 
is identically zero. Thus, we have obtained that 

EiWX 1 1|^ 0) _«] < cEpmo^X 1 )] < oo, 

which implies that ||X ||^o,-k < 00 almost surely, concluding the proof. 
Along the same lines as in the preceding steps, some I? bounds state that 

(78) E[\X\ S \\l s>ciBM ] < c(t - s) 2 «'o for < s < t < T 
and 

(79) E[\\Xi ||| S) £(B_ K ,B K )] < <t - s) 2 ^' for < s < t < T, 

for any k' q < H — 1/4 — k + v/2 and 7' < H, respectively. Following the 
same strategy as before, these bounds are enough to prove the remaining 
assertions of the proposition. □ 

Let us see now how this results can be related to our Hypothesis 1. Recall 
that the restriction k > 1/4 is dictated by the fact that we need to work in 
a space B K embedded in the space C([0, 1]) of continuous functions on [0, 1] 
in order to prove Proposition 4.10 below. 

Corollary 4.8. Suppose X is an infinite- dimensional fractional Brow- 
nian motion defined by the covariance function (43) for a given H > 1/2, 
with Q given by (44 ) for v > 0. Assume moreover that H > 7/8 — v/2. Then 
the incremental operator X 1 satisfies Hypothesis 1 for some 

k €(1/4,1/2), /t <H- 1/4- k + v/2, 7 = k + k. 

Proof. By the previous result, we have that X 1 has the required reg- 
ularity for any 1/4 < k < 1/2, kq < H — 1/4 — k + v/2 and 7 = kq + k < 
H — 1/4 + v/2. In order to check Hypothesis 1, we now need to require that 
7 + k > 1. In fact, there exists 1/4 < k < kq satisfying this inequality if and 
only if 7 + Ko > 1 , that is, 2H — 1/2 — k + v> 1. This is equivalent to assume 

H > 3/4 + k/2 - v/2 > 7/8 - v/2. 

In this latter case, it is easily seen that there exist 7, k, kq satisfying our 
requirements. □ 

Remark 4.9. If we are only interested in obtaining a local solution for 
our Young PDE, then the estimate (64) can be replaced by a bound in B K , 
which will be quadratic in y. Hence, using the fact that 

X 1 eCj£ K '~ K nC 2 K0 £ K ' K 
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for any 7 < H and 1/4<k<kq<H — 1/4 — k + P/2, the condition for the 
construction of the (local) fix-point map T becomes 7 + k > 1. To fulfill this 
requirement with our fractional Brownian noise, we only have to impose 
H > 3/4 — z//4. This condition is comparable (but a bit worse) with the 
results of [12], where the Hilbert spaces W 2k ' 2 were considered, and where we 
found H > 3/4 — P/2. One of the drawback of the approach presented in this 
paper is that the esimation of the random operators like X 1 in Banach spaces 
W 2a ' p for p > 2 seems very difficult. Moreover, it seems that the estimation 
in the Hilbert-Schmidt norm causes another small loss of regularity, which 
means that even in the case of a "regular" noise v = 1/2, our bound on H is 
H > 5/8 and not H > 1/2 as should be natural to expect and found in [12]. 
On the other hand, as we will see later, the operator approach seems better 
suited than the approach of [12] for a true rough-path expansion of SPDEs. 

Now that we have checked the assumptions on X , let us turn to the 
hypothesis on the nonlinear coefficient a in equation (69). In order to deal 
with the Sobolev norms, it is worth mentioning that, instead of working 
with the spaces B K = H K we have used so far, characterized by their Fourier 
decomposition, we will consider the Sobolev spaces W 2k ' 2 , induced by the 
norms 

(so) wm? - \Ml Hm) + / 0)1]2 

These spaces are obviously more convenient than the spaces B K for the 
computations on /, and they are closely related to these latter spaces, since 
the following classical relation holds true (see [1]): 

B K+£ C W 2k ' 2 C B K . e for any e > 0. 

Using these embeddings, we can consider the operator X£ s going from a 
space W 2k ' 2 to a space B K by just loosing a little regularity in t, s. Then we 
can verify that / satisfy a slight modification of Hypothesis 2. 

Proposition 4.10. Let a e C?(R) be a real-valued function. Then, for 
any k > 1/4, the function f:W 2K < 2 -> W 2k ' 2 defined by = <r(y(0) 

is locally Lipschitz, satisfies \f{x)\y/2 K ,2 < c/(l + |x|^2n,2) and is globally 
Lipschitz as a map f :B — >■ B. 

Proof. Recall that, for our particular situation, B = L 2 ([0, 1]), and it is 
easily checked that, whenever a G C 2 (M), the function / :B — > B is bounded 
and globally Lipschitz. 

With these considerations in mind, it is readily seen that / : W 2k ' 2 — >■ 
W 2k > 2 has linear growth. In order to check that / is also locally Lipschitz, 
note that its gradient can be computed as follows for y,h& W 2k ' 2 : 

V/(iO : W 2 ^ 2 W 2K > 2 , [V/(y) • h)(0 = a'(y(0)h(0- 
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Let us estimate now the norm (80) of V/(y) • h: first, if a' is a bounded 
function, then 

(81) ||V/(y) • /i|| L 2 ([0jl]) < ||er'|| 00 ||/i|| £ a([o ) i]) < ||o'|| 0O ||7i|| w *»,2. 

As far as the variational term of (80) is concerned, notice that we have 
assumed k > 1/4, which means that W 2k ' 2 C C([0, 1]), and for any h G W 2k ' 2 , 
WHoo < c\\h\\ W 2 K ,2. Thus, 

\[Vf(y)-h](0-lVf(y)-h}( v )\ 2 

[o,i] 2 ^ 



(82) 



1^(2/(0) -o'fofa))! 



2 



^ 1 1 t 1 1 2 ri i II "II 2 n ||2 1 

Putting together (81) and (82), we have thus shown that 

\\^f(y)\\c(w^) < + \\y\\ w ^), 
which easily yields that / : W 2k ' 2 — > W 2fi ' 2 is locally Lipschitz. □ 

Remark 4.11. Notice that, in spite of the fact that a is assumed to 
be a nicely behaved coefficient, its interpretation as an application from 
W 2k > 2 to W 2k ' 2 does not enjoy the usual assumptions of boundedness made 
on coefficients in rough path theory (see, e.g., [13, 16, 17]). This is one of 
the major sources of problems in our computations, and in general in the 
extension of rough path theory to SPDEs. 

Let us now summarize the considerations of the current section into the 
following theorem. 

Theorem 4.12. Let X be an infinite- dimensional fractional Brownian 
motion on [0, T] x [0,1], defined by the covariance function (43) and (44)i 
with H > 1/2 and v G [0, 1) such that H > 7/8 - 9/2 and let a G Cf (M). 
Then, there exists k G (1/4, 2H — 3/2 + v) such that for any initial condition 
tp G B K , the equation 

m = V(0. d t Y(t, = AY(t, dt + a(Y(t, Q)X(dt, d£), 

(83) 

t€ [0,TUG[0,1], 

with periodic boundary conditions, understood in the mild sense given by 
(57), has a unique global solution in C^' K . 
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Proof. By Proposition 4.10, the map / is Lipschitz and with linear 
growth from B K to £> K _ e for arbitrarily small s. As already noted this little 
mismatch of regularity can be compensated by the time-regularity of X . 
Then by a small modification of the arguments of Theorem 4.3 and by 
Proposition 4.4, we can directly solve the equation 



d t Y(t, = (A - ld)Y(t, dt + a(Y(t, Q)X(dt, d£), t€ [0, T] , £ € [0, 1] , 



as a rough evolution equation in C^' K . Now, if one wants to solve (83), it is 
sufficient to get an existence and uniqueness result for the equation 



which can be done along the same lines as for Theorem 4.3, by taking care 
of the additional drift term Ydt. This step is left to the reader. □ 

5. Rough evolution equations: the linear case. We pass now to the de- 
velopment of an expansion which allows to consider equation (54) in a case 
which goes beyond the Young theory, in terms of the Holder regularity of 
the driving noise x. We start with a simple linear case, that is, / = Id, which 
will hopefully lead to a better understanding of our method. 

5.1. Strategy. Recall that we wish to get some existence and uniqueness 
results for the equation 



Just like in the case of the Young integral, sketched at the beginning of 
Section 4.1, we will proceed as follows: 

(1) Expand (84) as if x were a regular process, until we get some terms 
which can be analyzed through the application of the operators 5 and 
A. 

(2) Define a natural extension of the notion of integral thanks to the first 
step, and show that this allows to integrate a reasonably wide class of 
functions. 

(3) Solve the equation in the sense given by this notion of integral. 

In the current section, we will mostly address the first of these three steps. 

If x is a regular process, equation (84) can be solved by means of the 
classical evolution theory. Furthermore, if y designates the unique solution 
to (84), then according to our expansion strategy, y also satisfies, for t, s G S2, 



3 t Y(t, £) = (A — Id)Y(t, dt + Y(t, dt + a(Y(t, £))X(dt, d£), 

*E [0,2*1, £e [0,1], 



(84) 
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However, the last term of this equation cannot be defined by applying the 
map A when x has low time regularity. In order to cope with this difficulty, 
let us expand again 8y by plugging relation (84) into the previous equation. 
Doing this twice, we get 



oyts — / Stu dx u S us y s + / Stu dx u I S uv dx v y v 

J s J s 

; pt pu 

Stu dx u S us y s -\- I Stu dx u / S uv dx v S vs y s 



I Sf u dx u I S uv dx v I S vw dx w S ws y s 

J s J s J s 

pt pu pv 

-\- I St u dx u I S uv dx v I S vw dx w 5y ws . 



Thus, going back to our notation on iterated integrals (41), we can recast 

(84) into 

(85) Sy = X 1 y + X 2 y + X 3 y + J(dx dx dx dx y), 

where, for i, s € £2 and 4> € B, the operators (^ l )j=i,2,3 are defined by 

(86) Xi<f> := Jt s {dxX'- l )cp = £ StudxuX 1 -^ 

with X® s = Sts- These operators are the new building block we will need in 
order to solve equation (84), and they play the role of the iterated integrals 
of rough path theory in our bilinear evolution context. Notice that the last 
term in equation (85) is considered as a remainder: suitable assumptions 
should be made to ensure that it will be small enough. Notice also that 
we stopped our expansion at the third order. We will see that this is the 
minimum order which allows to handle the Brownian case. 

Let us say a few words now about the algebraic properties of the operators 
X 1 : when x is a smooth process, we have, for example, 

pt pv pt pv 

^-^tUS / Stv dX V I Sy W dXyj S WS / Sfy dXy I S V yj dXyj Syj Ul 

J u J s J u J u 

and using some elementary algebra, we end up with 

pt pu pt pv 

^X-f- us / St V dXy S VU J S UU) dXyj Syj S ~\~ I S^ V dX v I S V yj dXyj S WU \_S US Id] 

J U J S J U Ju 

= Xtu-Xus + Xt u a us . 

Thus, taking into account our algebraic convention (16) and the definition 
of 5 given at (46), we have obtaind the relation 5X 2 = X X 1 . In a more 
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general way, it is not difficult to show by induction that 

5X n = Y J X t X n - i , 

i=l 

which are exactly the Chen relations in this setting. 

We can now specialize our previous program into the following: 

2a. Assume that the operator- valued 1-increments X ,X 2 ,X 3 are defined 
by some kind of operation which preserves the usual algebraic relations 
between integrals (e.g., use stochastic calculus with respect to an Hilbert 
space valued fractional Brownian motion or some other limiting proce- 
dure on discrete sums). They will be our (step-3) rough path. 

2b. Using (X , X 2 , X 3 ) define an integration theory for a sufficiently large 
class of functions Q so that it will be possible to give a meaning to 
integrals of the form z t = L St u dx u y u for any y € Q. We will call Q the 
space of paths controlled by X. 

3'. Study the map r : Q — > Q defined by T(y) t = L St u dx u y u , and prove 
that it has a fixed point y = T(y) which will be then a solution of the 
evolution problem (84). 

5.2. Integration of weakly controlled paths. We start by postulating some 
reasonable properties for X n . 

Hypothesis 3. We will assume that the process x allows to define some 
operator-valued increments X 1 , X 2 , X 3 , representing morally (49) and (86), 
respectively. This amounts for us to suppose that the X 1 's satisfy the alge- 
braic relations 

SX 1 = 0, 5X 2 = X 1 X 1 , 5X 3 = X 1 X 2 + X 2 X 1 , 

and that the following Holder-regularity properties holds true: 

x l e cl +{l - l)K0 c^'-f n d*°c™, * = 1,2, 3, 

for some r],p>0 and 7, kq such that 7 = ko + 77 + p and 7 + 3kq > 1. 

We will define now the class Q of processes we wish to be able to integrate 
against x: in the current situation, it will include any process which can be 
decomposed into a part depending on X l ,X 2 , plus a remainder term which 
is assumed to be small enough. For the sake of a contraction argument 
needed below (compare to the Young case), we fix a given time regularity k 
such that < k < kq. 
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Definition 5.1 (Weakly controlled paths). Let tp G B v be a given initial 
condition. A path y G C*' ri is said to be weakly controlled by X 1 , X 2 iiyo = ip 
and 5y can be decomposed into 

(87) 5y = X l y l + X 2 y 2 + y r , 5y l = X l y 2 + y 1 ^ 

with y % G C±' v i = 1,2, and a regular part y r G C^ 2 ^, y 1,r G C^ 4 "^ with 
k < kq Arj. Furthermore, we asssume that the regularity of y , y 2 and y r , y l,r 
can be related to those of X by the following relation: 7 + 3re > 1, a condition 
that can be always fullfilled by a suitable choice of k whenever 7 + 3ko > 1. 
Denote this space of controlled paths by Q K ,ri,ip > or when this does not lead 
to an ambiguous situation, simply by Q K)V or Q. Moreover, one can define 
a seminorm N on Q K-V in the following way: 

a% Qk] =N\yA' v ] + E -MM 00 '"] + E ^y l -A' v } 

i=l,2 i=l,2 

where we recall that the notation A/" has been introduced at Section 3.3. 

Remark 5.2. Even if a weakly controlled path is, strictly speaking, 
given by a tuple {y,y l ,y 2 ,y r , y 1,r ) we will, with a slight abuse of notation, 
denote it with its first component, that is, simply y. 

Remark 5.3. The notion of weakly controlled path appeared first in [6] 
in the finite dimensional context as a way to linearize the space of rough 
paths around the driving control. Even if this linearization does not preserve 
the whole structure of the space of rough paths, it is enough to find solutions 
of rough differential equations. 

With this notation at hand, we will try to implement now the strategy 
designed at the beginning of Section 4.1 in order to integrate a weakly con- 
trolled process y: let us first assume x is a smooth process, and y G Q. Then 
J{dxy) is well defined, and thanks to equations (48) and (87), we have 

J{dx y) = J{dx S)y + J(dx 5y) 

= J(dx S)y + J(dxX 1 y 1 ) + J(dxX 2 y 2 ) + J(dxy r ). 

Furthermore, for s < t, the term Jt s {dx X 1 y 1 ) above only involves y~, and 
hence the increment Jt s {dx X l y l ) is equal to Jt s {dxX l )y\, that is, X 2 s y\. 
This yields 

(88) J(dx y) = X x y + X 2 y 1 + X 3 y 2 + J(dx y r ). 
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Note that, in this last expression, the terms X 1 y, X 2 y 1 and X 3 y 2 are well 
defined under Hypothesis 3. In order to have a well-defined expression for 
J{dxy) in the rough case, it remains to handle the term J(dxy r ). Then let 
us write 

J(dxy r ) = J(dxy) - X l y - X 2 y l - X 3 y 2 , 
and let us analyze this relation by applying 5 to both sides. This gives 

(89) 6[J(dxy r )\ = -5[X l y] - *[*V] - o[X 3 y 2 }, 

and notice that in the last expression, 5\J(dxy r )] 7^0, since y r belongs to 
C-2 instead of C\. Moreover, a slight extension of Lemma 3.2 shows that, for 
M G d 2 (£(V)) and L G C X (V), we have 

8 (ML) = 5M L — M 5L = 5M L - M 8L. 

Applying this elementary relation to (89), we end up with 

S[J(dx y r )\ = -5X 1 y + X 1 ^- SX 2 y 1 + X 2 by 1 - 5X 3 y 2 + X 3 5y 2 

(90) = X\5y - X 1 y l - X 2 y 2 ) + X 2 ^ 1 - X l y 2 ) + X 3 by 2 
= X 1 y r + X 2 y 1 > r + X 3 5y 2 

under our hypothesis on y and X we have the following regularities: 
X l y r G Cj +3K '- p , X 2 y 1 ' r G Cj +K0+2 ^- p , X 3 6y 2 G C~< +2K0+ ^- p , 

so if 7 + 3k > 1 we can apply the operator A and express J(dxy) in terms 
of 5 and A only. Plugging (90) into (88), we get 

(91) J{dxy) = X l y + X 2 y 1 + X 3 y 2 + k(X l y r + X 2 y 1,r + X 3 5y 2 ). 

Similar to what we did in the Young case, we are now able to invert the 
procedure which lead to relation (91), by just invoking the assumptions 
made on X 1 and y: 

Theorem 5.4. Let x be a path such that X l ,i = 1,2,3, are well defined, 
and such that Hypothesis 3 holds true. Let also y G Qn,^,^ for < k < Ko < 
7 — k and K<n. Define z G Ci(B v ) such that zq = tp and 5z satisfies 

Sz = J{dx y) = X l y + X 2 y l + X 3 y 2 + k{X l y r + X 2 y 1,r + X 3 8y 2 ) 

and let z 1 =y, z 2 = y 1 , z l > r = X 2 y 2 + y r so that Sz 1 = X 1 z 2 + z 1 ' 1 ' . Then: 

(1) z is well defined as an element of Q Kri , and coincides with the usual 
Riemann convolution of y by x in case x and y are smooth processes. 
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(2) The seminorm of z in Q K ^ can be estimated as 

(92) M[z; Q K „] < c x T K °- K (\\ ^|| B „ + M[y; Q KiV }), 

for a positive constant cx depending only on X % , i = 1,2,3. 

(3) It holds that, for any < s < t <T 

n 

(93) J ts (dxy) = lim V S tu+ M+uUVU + K+^A + X l +1 ,tA ], 

|n ts |->o . =Q 

where the limit is over all partitions Ut s = {to = t, . . . ,t n = s} of [s, t] as 
the mesh of the partition goes to zero. 

Proof. We will divide again this proof in several steps. 
Step 1 : Let us start by evaluating the regularity of the terms in the right- 
hand side of (91), that is, 

A = X 1 y, B = X 2 y\ C = X 3 y 2 , 
D = A(X 1 y r + X 2 y 1 > r + X 3 Sy 2 ), 

under our standing assumptions. 

In order to bound A, we will first estimate |y s |e itself for s <T: if y G Q v , 
we have y s = S s ip + 5y s o, and hence 

(94) \\y\\o,B v <\MB v +T K M[y;C^}. 

In particular, y is bounded in B v on [0, T]. Thus, if X 1 G nC^C™, 
we have X 1 y G Cj' f \ an< ^ a ^ so X V G C^ ' 71 ■ Moreover, 

\Xly s \ Bv < \\X l \\ KQ ^(t-sY°(\^\ Bn +T*N[y^}), 

and thus 

(95) NiX^C^] < \\X l \\ K0 ,vM\B v + T^[y;Cl' v ])T^. 

Let us estimate now the term B, that is Af[X 2 y 1 ; since y 1 G C^ ' 11 

and X 2 G C 2 " ^, we obtain again that XV G (^ko+k,^ and we have 

(96) M[X 2 y y -Cl^\ < l^h^Mly 1 ^]^-^. 



The term C can now be bounded along the same lines as for A and B. 

3 

assumed that 7 + 3k > 1. Thus, the operator A can be applied to X y r + 



Moreover, for the term D, as we already observed above, X l y r € Cg -1-3 '"' p , 
X 2 y^ r G c7 +K0+2K '- p and xHy 2 G c7 +2K0+K '- p , and observe that we have 
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X 2 y 1 ' r + X 3 8y 2 , and invoking inequality (40), we get that 

\\k{x 1 f + x 2 y^ + x 3 5y 2 )h+^- P 

< c\\x l f + xV' r + x 3 5y 2 h + ^- P 

(97) 

< c(\\X%^ p N[y r A K "\ + UX+^-pM^'A^] 

+ \\X% +2Ro ^ p M[y 2 ;C^\)- 

Summarizing inequalities (94)-(97), we have obtained that z is a well-defined 
element of C^' v , and that it satisfies 

\\5z\\ K , v < c x T K °- K m Bv +M[y, Q K , V ]). 

Step 2: Let us estimate now z as an element of Q K ^ V . The natural decom- 
position of 5 z is obviously Sz = Xz 1 + X 2 z 2 + z r , with 

z 1 =y 1 z 2 =y 1 and z r = X 3 y 2 + k{X l y r + X 2 y^ r + X 3 5y 2 ). 

It is now easily checked, along the same lines as for Step 1, that z satisfies 
relation (92). 

Step 3: In order to see how to get the convergence of the Riemann sums 
to J{dxy) it is enough to remark that 5z can be written as Sz = (Id — 
A5)[X 1 y + X 2 y 1 + X 3 y 2 ]. Applying Corollary 3.6, we now get relation (93). 
□ 

Remark 5.5. The space of weakly controlled paths is a vector space 
with respect to the action of M but not with respect to other interesting linear 
endomorphisms of B. The problem lies in the fact that for general linear 
L:B^B we can have 5Ly ^ L5y since L does not necessarily commute 
with the semigroup (which appears in the definition of 5 = 5 — a). 

5.3. Linear evolution problem. Let us turn now to the main aim of this 
section, which is to get an existence and uniqueness result for equation (84). 

Theorem 5.6. Assume that Hypothesis 3 holds for the triple of in- 
cremental operators X X ,X 2 ,X 3 with 7, kq,k,tj,p such that 7 = kq + 77 + p, 
7 + 3kq > 1 and K < kq . Then: 

(1) Equation (84) admits a unique solution y G Q v . 

(2) The map (ip , X 1 , X 2 , X 3 ) ^ y is continuous. 

(3) For (t, s) £ S2, the map & t s ■ B v — > B v , such that Qtstp = Vt when y s = ip 
and 5yts = Jt s {dxy) is a bounded linear endomorphism of B v , and it 
satisfies the cocycle property & t u&us = &ts- 
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Proof. Like in the Young case, the solution y will be identified as the 
fixed point of the map V : Q Kjr) — > Q Kjr) defined by z = T(y), with zq = ip and 
fizts = Jts{dxy). And here again, we will concentrate on the fact that, on a 
small interval [0, T], the ball 

B = {y,y = 4>,N[y,QK,r 1 ] < h%J 

is left invariant by the map F. 

Indeed, whenever y € B, then Theorem 5.4 asserts that for z = T(y), the 
following estimate holds true: 

M[z; Q K>V ] < c x T^ K (\^\ Bri +M[y; Q K , V ]). 

Hence, if one chooses a small enough T, so that cxT K °~ K < 1/2, it is readily 
checked that M[z; Q« )7? ] < which proves that z € B. The contraction 

property is now a matter of standard arguments, and the remainder of the 
theorem follows easily. □ 

5.4. Application: stochastic heat equation. In the sequel of the paper, 
for sake of simplicity, the generic situation of a process X with 7-Holder 
continuity in time with 7 < 1/2 will be the case of an infinite-dimensional 
Brownian motion, given by the covariance function (45). For this special 
process, we will try to construct a pathwise solution to the linear stochastic 
heat equation on [0,1]. At the end of the section, we will give some hints 
about the way the fractional Brownian case should be treated. 

5.4.1. The Brownian case. Like in the Young case, the key step in order 
to apply Theorem 5.6 to the Brownian setting is to define (^*)i=i 2,3 m a 
reasonable way, and then to check Hypothesis 3. We have chosen here to 
deal with an ltd type definition for X n , and we get the following result: 

Proposition 5.7. Let X be an infinite- dimensional Brownian motion 
defined by the covariance structure (45), with Q given by (44) for v £ [0,1]. 
For n = 1,2,3, let X n be the incremental operators given by (70) and (86), 
respectively, where the stochastic integral has to be understood in the ltd 
sense (see, e.g., [3, 25] for a complete definition). Then, almost surely, 

v-n c /?7+(n-l)«o r r),-p n Ann r t],T] 
fc u 2 ^HS 1 1 u 2 ^-HS 

for any n > 1/4, 7 > kq> k satisfying 

k < 1/4 — 77 + P/2 and 7<l/2, 

with v = inf (y\ 1 /2) . 
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Proof. We have already proved the regularity of X in the fractional 
Brownian case. The proof in the current case would be similar, and we omit 
it. It will be enough to take 77 = 1/4 + s, kq = u/2 — 2e and 7 = kq + p + e, 
p = 1/4 — v/2 for a given small e > 0. 

Let us concentrate then on the regularity properties of X 2 : we will prove 
in fact first that X 2 € CT^^ C^a p , and for this, we will proceed along the 
same lines as for the proof of Proposition 4.4. Let us sketch the main steps 
which have to be followed. 

Step 1: First of all, we have to estimate \\Ao p X 2 s A~ v \\iis-i3^ > .b, and it is 
readily checked that Ao P X 2 s A~ v is represented by the kernel 



Thus, when considered as an operator from to B- p , we obtain that 



J[0A] 2 

Moreover, some standard considerations about iterated integrals for Brow- 
nian noises (see, e.g., [3, 25]) yield 





(98) 





with 




Plugging this equality into (98), we end up with 



£[II^ 2 s |Ihs]= / du / G-^Jm^ui^m) 




x Q(m -fji)dr]idfji, 



where we have set e = t — s and 




x G 2v v (r] 2 ,fj 2 )Q{m ~ m) dr] 2 di) 2 . 
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Furthermore, using the spectral decomposition of Gt and Q, introduced 
respectively, by (42) and (44), we obtain 

G u - v (rn,r]2)G u - v (f)i,f)2) G^ v (m , m ) Q (V2 - m ) drj2 dfj 2 

[o,i] 2 

e -{^i+^i)[u-v] e -2\ k v 

= 2^ e i(jll) e j(f)l) -—2^ 1 {i-k-l=0} 1 {j+k+l=0}- 

i,j,k,leZ I k 

Injecting again this value into (99) and using the fact that \-i = A.;, we have 
that 

1 



^[II^sIIhs] - y, 



i,j,k,l,m,n£E A l A ™ A n A l A k 



due- 2Xm(£ - u) / dve- 2X >- v] e- 2XkV , 



with 

(100) E = {j,k,l,m,n eZ;m + n = j,k + 1 = —j}- 
Thus, we get 

£[H*t 2 s |lHs]= V ■> 1 9 [ £ du f 



1:^ 



1 F , P , , dv 

<c > — — 7T~ \ n du dvdu- 



'^neE^^K^KJo Jo ' (e-uHu-vf 
< ce 2 ~ a ~ 6 , „ t o du dvdu- -r. 

j,k,l,m,ndE J "* ra fc i 

The double integral above is finite whenever a, 6 € (0, 1), while the sum can 
be handled along the following lines: first, rewrite 

9= V 1 V 1 

m,n,j:m+n-j=0 A ra k,l:k+l=-j k A Z 

and observe that, thanks to Lemma 4.6 and according to our hypothesis 
rj > 1/4, we have 

k,l:k+l=-j A k A i k,i.k+l=-j A k A l 

where C stands again for a positive constant which can change from line to 
line. Then 

s<c y j-- — <cy^— y -JL_ 

m,n,j:m+n-j=0 A j Am A n m Am n,j:n-j=-m j « 
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and choose b > 1/2 — 2v, so that another application of Lemma 4.6 gives 

S<C^l/A^ +a+6+2p - 1 / 2 . 



rn 



This latter sum is finite when a + b> 1 — 2p — 2v. Then for any 9 <9* such 
that 29* < 1 + 2p + 2v, we have found that E[\\X%, ||| s ] < ce 2 -"- 6 < ce 2e . 

Step 2: One can go from L 2 to L p estimates for m just like in Proposition 
4.4 Step 4: indeed, we have 

■V]* f = i 

Moreover, Kt s (^i,rii) is a variable of the second chaos H2 with respect to the 
Gaussian field X, and invoking [19, Relation (1.61)], the L 2 and L 2p norms 
on H2 are equivalent. Thus, for any integ 6r p ^ lj there exists ci constant c^, 
such that 

E[\\Xl\\%]<c p {t-s)^. 
Step 3: We will conclude now thanks to Lemma 3.8, which reads here as: 

Il^ 2 |l72,»?,-P — c [U~/2+2/p,P,V-p(X 2 ) + Il^ 2 |l72,»?,-p] 

= c[L r 72+2 /p,p,r ) ,-p(^ 2 ) + H^ 1 ^ 1 H-ya^-p], 

for any integer p > 1. According to the previous step, it is then easily 
checked that, for any 72 < 9* = 1/2 + v + p, and p large enough, the term 
U~/ 2 +2/p,p,r),~p(X 2 ) can be bounded almost surely by a finite constant. Recall 
now that we have chosen 7 = ko + p + 77, /-to = — 2e and 77 = 1/4 + e. 
Thus, 

7 + k = 2k + p + r/ = P - r/ + 1/4 + p < 6>* , 

and hence U li+ 2/ P) p )rj - p {X 2 ) < 00 for any 73 < 7 + kq. 

Let us treat now the term XX. Along the same lines as in Proposi- 
tion 4.4, it can be shown that X 1 G & 2 CP >~p and X 1 G C%°C™. Hence, by 
composition of operators, we get X 1 X l G C^ +K ° C^" 9 , which means that 
IIA^X 1 || 72fli _ p is finite for any 73 < 7 + kq. Summing up this short discus- 
sion, we have obtained that 

||X 2 || 72jr)i _ p finite a.s. for any 73 < 7 + kq. 

One can proceed then to prove that X 2 gC 2ko />'' ? by a slight elaboration 
of the computations above. This easy exercise is left to the reader. 

The proof for the operator X 3 follows the same lines and will not be 
reported. Indeed, we prefer to concentrate on the regularity properties of 
higher order operators in the more complex situation of Section 6.4. □ 
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We are now able to apply our abstract results to the stochastic heat 
equation. 

Theorem 5.8. Let X be an infinite- dimensional Brownian motion on 
[0, T] x [0, 1], defined by the covariance function given by (45) and (44) with 
v > 1/3. Then there exists 77 > 1/4, 0<k<7<1/2 such that k < kq and 
7 + 3k > 1 such that, for any tp E B v the equation 

Y(0, = V>(£), d t Y(t, = AY(t, dt + Y(t, QX(dt, d£), 

ie[0,T],£e[o,l], 

with periodic boundary conditions, understood as equation (84), has a unique 
solution in Q K ^,^- 

Proof. Like in the proof of Theorem 4.12, the claim is readily checked 
once we have shown that X n , n = 1,2,3 satisfy Hypothesis 3. This amount 
to check that there exist kq > k and 7 < 1/2 such that 

k < l/4-77 + z?/2, 7 + 3k>1. 

However, thanks to Proposition 5.7, it is enough to take r\ = 1/4 + e, kq = 
v /2 — 2e, k = kq — £, p = 1/4 — v/2 and 7 = kq + p + ij for some small e > 0. 
The condition 7 + 3ko > 1 can then be read 7 + 3ko = 1/2 + 3P/2 — 4e > 1, 
which is possible whenever v > 1/3. □ 

5.4.2. The fractional Brownian case. In order to define an integration 
theory for the fractional Brownian motion beyond the Young case one has 
to start, like in the Brownian case, by defining the operators X 1 and X 2 in 
a natural way. We have already seen that X 1 could be understood by means 
of Wiener integrals, and for X 2 , two reasonable choices for the definition of 
(86) seem to be the use of either Skorokhod or Stratonovich integrals with 
respect to the fractional Brownian motion X. However, it turns out that 
these two solutions are equally unsatisfactory, for two different reasons that 
we proceed to detail now: 

(1) When one computes moments of random variables of the second chaos 
defined by Stratonovich integrals, some trace terms appear, a classical 
phenomenon which is explained for instance in [19] in the general case, 
in [23] for the stochastic heat equation, or in [21] for the fBm. In the 
current situation, if we want these trace terms to be convergent for a 
fractional Brownian motion X defined by (43) and (44), one has to 
choose v > 1/2, which means in particular that Q is a bounded function 
of £ G [0,1]. In other words, we are not allowed, even if H > 1/2, to 
consider a distribution-valued noise in space, which was one of our main 
aim. 
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(2) The Skorokhod integral works better as far as convergences and regu- 
larity estimates are concerned. But one of the basic ingredients of our 
algebraic manipulations on integrals is the fact that one can write, under 
suitable hypothesis: 



Stu dx 



s 



S tu dx 



an equality which is known to fail in the Skorokhod case (see [19] again 
for further explanations). For instance, the relation 5X 2 = X 1 X 1 , which 
is useful in our analysis, does not hold true when 5X 2 is understood as 
a Skorokhod integral. 

In order to cope with these problems, one can adopt the following strategy: 
compute the correction term P, understood as a 2-increment operator- valued 
process, which allows to write 

(101) 5X 2 = X 1 X 1 + P, 

when X 2 is defined via Skorokhod integration. Notice that, since X 2 is an 
element of H2, the process P is deterministic. 

Recall now that the operator 5 has been defined as follows: for a suitable 
Banach space V and M £ C9°{V\ set 

(102) 5M = 8M - Ma = 5M - aM - Ma. 

Then the operator 5 enjoys the same kind of properties as 5, and in par- 
ticular, 55 = and ker^ = Im5\^. Moreover, relation (101) can be read 

as 5P = 0. Thus, there exists another process T € C2C such that ST = P. 
Consider then an explicit version of T and set X 2 = X 2 — T. Then X 2 is still 
a Levy area type process, such that 5X 2 = X 1 X 1 , which means that hope- 
fully, X 2 will enjoy both algebraic and analytic properties allowing a nice 
extension of the notion of convolution integral. However, an open problem is 
to understand in which sense the integrals defined using this corrected Levy 
area X 2 are useful and/or natural. We plan to report on this possibility in 
a further paper. 

5.5. The algebra of a rough path. Bilinear stochastic equations, in finite 
or infinite dimensions, are often handled by means of chaos decomposition 
(see, e.g., [14, 15]). In this section, we will try to stress some relationships 
between our pathwise approach and this latter method. 

Our basic Hypothesis 3 states that 

71-1 

(103) 5X n = Y J X k X n - k 

k=l 
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for n = 1, 2, 3, and moreover that 

(104) X n e c<y+(n-l)K0 £ r,,-p and X n £ C™° , 

with 7 + 3ko > 1, 77 > 1/4 and p = 1/4 — P/2. Furthermore, it can be shown, 
along the same lines as for Theorem 3.5, that there exists an inverse A to 8 
on C^C°~ P fl ker(5) for a certain /x > 1. 

Let us see now how to construct an operator X 4 satisfying the operator 
Chen relation (103): by composition of operators, and invoking Hypothe- 
sis (104), it is easily checked that X 1 X 3 + X*X l + X 2 X 2 G C^ +3k ° C r > <~p. 
Furthermore, we have assumed that 7 + 3kq > 1, and thus, by analogy with 
Theorem 5.4, we will set now X 4 := AfA^X 3 + X 2 X 2 + X 3 X 1 ], which is 
well defined as an element of S]^ /^ >~p and thus that belongs to C^C™ 
(since 7 = kq + 77 + p). It turns out that this procedure can be iterated, and 
we obtain the following proposition. 

Proposition 5.9. Let X satisfying Hypothesis 3. Then one can con- 
struct a sequence {X n ;n > 4} out of X l ,X 2 ,X 3 , such that, for any k < 
Ko = ln - P} W e have X n e s] +(n ^ 1)K °a ~P, 

\\x n \\cr°c^p^ c ( nl )~ K0 

and such that the operator Chen relations (103) are satisfied. 

Proof. The X n are constructed by an induction on n. Then it is clear 
that X n G £;7+(«-i)«o £ r,,- P _ Moreover, for n > 4 we have ukq > 1, so that 
X n = A(^fc= 1 1 X n ~ k X k ) can be defined directly as an element of C™ . Then 
the same kind of arguments as in the finite dimensional case [10] prove that 
we have the inductive bound 

(105) \\X n \\ d n KOc ^<C x (nl)- K0 . □ 

In such a setting, the lifted rough path allows to express the ltd map 
which sends initial conditions to solution of the linear equation (84) by a 
convergent series of operators. 

Corollary 5.10. Under the conditions of Proposition 5.9, there ex- 
ists an operator T , defined as an element of £^ +(n ' 1)K ° C^'~p , given by the 
strongly convergent series T := ^2^LiX k , and such that the solution of the 
linear problem (84) satisfy the equation Sy = Ty, or written in another way 
Vt = S ts y s + T ts y s , (t,s)eS 2 . 

In particular, if we define Xf s = St s and set T = X° + T we have that T is 
a cocycle of operators: 

T t s = T tu T us , (t,u,s)eS 3 . 
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Proof. The convergence of the series for T in the operator norm follows 
from the bound (105) on ||X n ||. The cocycle property is proven as in finite 
dimension. The uniqueness of the solution to the linear problem allows to 
identify the operator T as the ltd map for the rough evolution equation. □ 

6. Polynomial nonlinearities. Going back to the general setting explained 
at Sections 3.2 and 3.4, we will consider now an equation of the form 

(106) Vt = Stil>+ f S tu dx u M n {yT), 

Jo 

where y lives in the Hilbert space £>, where M n :B n —±Bis some unbounded 
multilinear operator from the Hilbert tensor B n = B® n to B, and where we 
understand <p® n = ip ® • ■ • <8> (p € B n as the tensor monomial generated by 
ip € B. In fact, for sake of simplicity we assume that M n is symmetric and 
we restrict our discussion to the case n = 2, letting M<i = B, the general 
situation posing no more conceptual difficulties. Then our general strategy, 
like in the linear case, will be first to expand equation (106) for a smooth 
driving process x in order to guess the appropriate rough-path underlying 
this equation. It will be seen that those expansions involve some increments 
indexed by trees. Studying the algebraic and analytic properties of these 
increments, we shall obtain a reasonable notion of solution to our quadratic 
equation. 

6.1. Formal expansions and trees. Let us first simplify a little our setting. 
Recall that we wish to solve an equation of the form 

(107) Vt = Stfl>+ /"* S tu dx u B(yf), 

Jo 

where we specialize our situation in the following way: assume first that B = 
L 2 ([0, 1]), which means that we are back again to the heat equation setting 
of Section 3.4. Then B : B® 2 -> B is defined by [B{4> ® VOKO = HO^PiO for 
£ £ [0, 1] , whenever this expression makes sense in B. Assume for the moment 
that x G C\C K,K for k large enough. We can expand equation (107) as 

yt = S ts y s + S tu dx u B((S us y s )® 2 ) 

J s 

+ 2^ S tu dx u B^S us y s ® j S uv dx v B((S vs y s )® 2 )^ 
+ J Studx u B(J S uv dx v B{{S vs y s )® 2 ) 

S U v dx v B((S vs y s )® ) 
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(108) 

+ 4 / S tu dx u 



x B \ySusVs ® J S uv dx v 

x B \{S vs y s ) £ S vw dx w B((S ws y s )® 2 )^ \ 

+ h.o. iterated integrals. 

As we see, iterated integrals appear here in combinations which are not as 
easy to handle as in the bilinear case of Section 5.3. A natural way to code 
this kind of expansion is to use planar trees, as explained below. 

Without entering too much into formal definitions involving trees, let us 
mention that we shall consider planar binary rooted trees T of the form 



v 



X/,xy, v ,y,\/,^>,y,etc. 



which allow to give a compact expression of the iterated integrals appearing 
in the expansion (108). Observe that each tree can be constructed from the 
trivial tree To = • by using the binary operation V : T x T — > T consisting 
in gluing two trees at a newly created root, so for example 

•••• 

V V 

V =V(V(t ,t ),V(t ,t )). 

For any tree r £ T, we associate the function d(r) that counts the number 
of leaves on the trees, so that cZ(to) = 1 and d(V(ri, T2)) = d{r{) + dfo). 
Let us see now how to represent expansion (108) thanks to planar trees. 

Define recursively an operator- valued increment X T € ^CiB^^ ; B) for r € 
T as 

ft 



(109) X£=S ts and xj}^ = ! S tu dx u B(X T u \® X 



T2 ) 



Notice that X T has always to be considered as an operator acting on B^ T \ 
For instance, we understand that, if r = V(tx,T2), we have 

xl^\ Vl ®---® Vd{T) ) 

i-t 

Studx u B(Xl\((pi ® ••• ® <Pd( n )) 

This latter formula justifies also in a sense the use of planar trees, since in 
general X V ^ T1 ' T2 ' > X V ( T2 ' T1 \ In order to illustrate this fact, consider the 
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simple example where t\ = tq and t 2 = V(tq,tq). Then d(V(Ti,T2)) = 3 and 
^ (T1 ' T2 V ®W® Vs) = / S tu dx u B(Xll(vi) ® xl^ ro \^ 2 ® 

J s 

while 

^ (T2 ' T1 V ^^2 0^3) = f S ta dx v B(X$™Xvi®<p 2 )®X£(<p 3 )), 

J s 

which are a priori clearly different objects. 

With this notation in hand, it is now checked that our previous expansion 
(108) can be written in a simpler way as 

(110) (Sy)ts = X? s (yf 2 ) + 2XZ(yf 3 ) + X^(yf*) + AX%(yf 4 ) + r, 

where r € C2C8) is some remainder term, and where we took care to dis- 
tinguish the various operators obtained by permuting the factors in the B- 
tensors. Of course, we could have expanded the solution further, and some 
operators associated to larger trees would have appeared. However, in a 
smooth enough situation, the strategy in order to solve (107) is now clear: 
we can use the map A to eliminate the remainder from the equation: 

(111) 8y = (1 - k5)[X~(y m ) + 2X^(y m ) + X^(y m ) + 4X^(y® 4 )] 

and try to solve this by fixed-point method. The only condition we need to 
check is that 

(112) S[X"(y® 2 ) + 2X^(y® 3 ) + X v - V (y 04 ) + AX^{y m )] 

should be in the domain of A, which means that its time-regularity should 
be greater than 1. The computation of expressions like (112) requires a little 
algebraic preparation. 

6.2. Algebraic computations. To ease some computations, we introduce 
an "improper" increment Ets = Id (improper because it does not vanish as 
t = s), where the Id has to be understood, according to the context, as the 
identity operator on the vector space under consideration. For example, we 
can write 5h = 5h — (S — E)h. Moreover, we also introduce et s = 1 taking 
values in R, so that for example, if z G then (ze)t us = zt u e U s = z tu . 

It will also be useful to extend the action of S to the tensors B n by letting 
5z = 5z — (S — E)z, where S :B n — > B n is defined as S = S <g> • • • ® S for any 
n > 1. If the reader is unconfortable with giving the same name at differ- 
ent operators, he can think that S is defined on the direct sum © n>1 <B n ; 
furthermore, we will write explicitly S n when the context is insufficient to 
determine the actual space on which S is defined. Analogously to the case 
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of 8, the operator 8 defined by (102) can be allowed to act on CiC(B n ,B) as 
8H = 8H - (S - E)H - H(S - E). 

We wish first to understand how the operators 8 and 8 act on tensor prod- 
ucts. More specifically, we shall need three relations which are summarized 
in the following lemmas. 

Lemma 6.1. The following relations hold true: 

(1) Let z,w GCi(B). Then 

(113) 5(z®w) = Sz®8w + 8z® Sw + 8z®8w. 

(2) Let z,w eC 2 (B). Then 

8(z®w) = ze®8w + ze® Sw + Sz ® we + Sz ®we 
(H4) . „ 

+ Sz (g> 8w + 5z <g) Sw + 8z® Sw. 

(3) Let Z G C 2 C(B a , B b ) and W G C 2 £(# c , £ d ). T/ien 

8(Z ®W) = ZS a ®~8W + ZS a <g> S d W + 8Z® WS C 

(115) . . . „ 

+ S b Z ® WS C + S b Z ®5W + 8Z® S d W + 8Z® SW, 

where an example of notational convention is given by 

(ZS a ® s d w) tus = (z tu sZ a ) ® G £(S a+c , £ c+d ). 

Proof. These relations are easily checked by elementary computations. 
We include the proof of the third one for sake of completeness: notice that 

SZ = Z° — ZE - EZ — (S b - E)Z - Z(S a — E) = Z° — ZS a - S b Z, 

where Z^ us = Z ts , and thus 

8(Z <g> W) = Z° <g> W° - S b Z <g> S d W - ZS a ® WS C 

= (S b Z + ZS a + SZ) (S d W + WS C + SW) 

- s b z o s d w - zs a ® ws c , 

which yields relation (115) by a straightforward expansion. □ 

We also want to understand how S,S act on the operators X T . A first 
relation in this direction is to note that, according to Lemma 3.2, if X T G 
C 2 C{B® n ; B) and h G C^B® 71 ), 

(116) S[X T h] = (SX T )h-X T 5h, that is, S[X T h] = (SX T )h - X T Sh. 



58 



M. GUBINELLI AND S. TINDEL 



It is thus useful to compute quantities of the form 5X T . To this purpose, 
consider n > 1 and define / : C 2 C(B n , B 2 ) -> C 2 C(B n , B) by 

(117) I(H) ts = f S tu dx u BH us . 

J s 

This kind of expression can be related to our tree-indexed increments by 
noticing, for instance, that 

X~ = I(S 2 ), X* = I(X"®S 1 ), X* = I(S 1 ®X~), 

and generally speaking, (109) can be read as 

x V{r U T2) = j^n ® X T2 ), X v{ - T °'^ = I(S ® X T ) and 

(118) 

X nr,r ) =I (X T ®S). 
Hence, we shall compute differentials of terms of the form 1(H). 

Lemma 6.2. Let H <^C\L(B n ,B 2 ). The following formulae hold true for 
the derivative of 1(H): 

(5I(H)) tus = L tu (S 2 )H us + I tu (5(H)) 

and 

(5L(H)) tus = Itu(S 2 )H us + I tu (S(H)). 
Furthermore, if we assume that 5H can be decomposed as (5H)t us = 
Y,j<M H iu j)H us j \ for a given M > 1, H^ e6\C(B 2 ,B 2 ), and H^) e 
C\C(B n ,B 2 ), then we obtain 

51(H) = I(S 2 )H+ ^ I(H^)H^\ 

Proof. We have 

[5I(H)] tus = Its(H) — Itu(H)S us — StuIus(H) 

t rt 

&tw dx w BH WS / St w dx w B H WU S US 



Stu I S un , dx w BH i 



tu J >~>uw u x w u ±i ws 
t rt 

&tw dxyj BH WS / St w dx w B H WU S US 
J u 

t 

Stw dx w B(S WU H US ) 
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+ / Stw dx w [BH WS — B(S WU H US ) — B H WU S US ] 

J u 

= I tu (S 2 )H us + I tu (5(H)), 

which proves our first assertion. The second one is now trivially deduced. 

□ 

With these preliminaries in hand, we can now compute the action of 5 on 
the tree-indexed increments we have met so far, in the following way. 

Lemma 6.3. Let x be a smooth C(B) -valued path. Then we have 
5X^ = 0, 

5X^ = L(S 2 )(X~ <g> S) = ® S), 

6X"- V = X V (X V ® X v ) + X*(X~ ®S) + X^(S ® X v ) 

and 

5X^ = X V (X* ®S) + X*(X" ® S 2 ), 

5X* = X^{X* ®S) + X^(S 2 ® X v ), 

5X* = I v (S®I x ')+I ,y (I v ®5 2 ), 
5X^ = X~(S X*) + X' y (S 2 <g> 

Proof. All these relations are obtained by elementary computations, 
and we shall only sketch the proof for some of them: first of all, invoking 
Lemma 6.2, we get: 

5X^ = 61 (S 2 ) = I{S 2 )S 2 - I(S 2 )S 2 = 0, 

where we used the fact that 5S = —SS. As far as X s * is concerned, we have 

6X* = 5I(X- ® S) = I{S 2 ){X- ® S), 

since <5(X V ® S) = by a direct computation using formula (115). Similarly, 
it holds that 

5X V - V = 5I(X~ ® X") 

= /(<S 2 )(X V ® X^) + I(S ® A V )(A V ®S) + I{X^ ®S){S® X^), 

owing to the fact that 

5(X~ <g> A v ) = 5X V <g> X~S + X~S ® 5X V . 
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Now, invoking (118), we have I(S <g) X^){X^ ® S) = X* and I(X^ ® S) = 
X^, which yields 

SX v - y = ® X v ) + X^(X V ®S) + X^(S <g> X v ), 

as claimed in our lemma. □ 

It is important to note now that all the previous computations have been 
performed for a smooth path x. However, we shall ask our driving process 
x to satisfy the following assumption: 

Hypothesis 4. We assume that the path x allows to define some in- 
crements X T for any r € T such that d{r) < 4. We also suppose that those 
increments satisfy the relations of Lemma 6.3, and that the following Holder 
regularities hold true: setting \r\ = d(r) — 1, we have X T € C^ K ° C(B^ T \ B^) 
and X T eC^ +{ ^-^ K0 C(B^ T) ,B- P ), withj + nK >l and 7 = n + n + p, for 
a given n > 1/4. 

Remark 6.4. Here again, it is important to work in spaces of the form 
B v with n > 1/4. Indeed, these spaces are algebras, which ensures at least 
that, whenever (f>,ip & B v , then B((f>, ip) € B. 

Remark 6.5. The peculiar relation between the various parameters in- 
volved in Hypothesis 4 has been suggested by the example which we treat 
later on and it is due to the mixing between space and time regularity due 
to the analytic semigroup. As operators the X T, s can map to more regular 
spaces (with respect to the scale associated the generator A) at the price 
of loosing some time regularity. This is a phenomenon which is not pecu- 
liar of infinite-dimensional rough paths associated to random processes but 
it is found also in the rough-path approach to deterministic PDEs like the 
Korteweg-de-Vries equation or the Navier-Stokes equation [7-9]. 

6.3. A space of integrable paths. The general discussion of the bilinear 
equation requires a deep understanding of the algebra of X. It is not the 
aim of this paper to enter into these kind of considerations, and we prefer 
here to concentrate on a particular case where k is sufficiently large to stop 
the expansions at some low (but nontrivial) order. So here we assume that 
7 + 3k> 1. 

In order to solve the fixed-point problem associated to (107), we introduce 
a new space of weakly controlled paths, denoted by Qx, K , which enjoys some 
nice stability properties under the map r : y 1— > z = Sip + I(y <S> y)- 
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Definition 6.6. Let tp G B^ an initial condition, and x a driving noise 
satisfying Hypothesis 4, with 7 + 3k > 1. We say that a path y G Cl' K (B ri ) 
belongs to Qx,k if 2/o = V'j an d can be decomposed into 

(119) iy = IY + ^V + ^Y + yK 

where y^y* and y^ can be written as: 

y v = w ®w, y~^ = w^<S>w, y <y = w<giw s ', 

and the following regularities hold true: 

y G Cl' K (B v ), w,ye C^(B V ), G 6*^(8$), y» G Cf 2 (^), 

where k > k\ > K2, k — ki = k\ — K2 = n and 7 + 3/^2 > 1- On <2x,k> we define 
the seminorm 

A% Q x ,«] = A%C 1 " (£„)] +N[w-C^{B V )} 

+ Af Cf 2 (Bj)] + N[y* ; C 2 3 « 2 (B,)] . 

Note that the constant path yt = Stip is a controlled path whenever ip G 
B K , and in this case w,w v ,y^ are all identically zero. Furthermore, the space 
Q satisfies the following useful stability property. 

Theorem 6.7. Assume that x satisfies Hypothesis 4> where we recall 
that 7 + 3k > 1 and kq > k. For y G Q K ,x, define z = r(y) G Q K ,x by zq = ip 
and a decomposition of the form 

5z = X v z v + X^z^ + X^z* + z*, 

with z w = w z ® w z , z^ = w z £>5>w z , z* = w z ® w^, 

where w z = y, = y v , and z$ G C^'iB) is a remainder which can be written 
as 

z« = X v - V (y v <g> y v ) + X^(y^ <g> y) + X^iy* ® y) 
+ X*(y ® y*) + X^(y ® y*) - A( J), 
where J is defined by relation (121). Then: 

(1) T : Q Kj x -> Q Ko ,x is we// defined. 

(2) coincides with I(y (g> y) in t/ie smooth case. 

(3) TTie following estimate holds true: for all < S <T we have: 

M[z;Q x , Ko ([0,S])] 

(120) 

< c x (i + ivk + l^ok + H\v + §MA % Qx,k([^s})})\ 

for a positive constant Cx which only depends on the rough path X . 
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Proof. Start from two smooth paths x and y. If we apply the I map 
defined at (117), we obtain, just as in the Young case (48), 

5z = I(y ® y) = I{6(y ® y) + S 2 {y ® y)) = I{S 2 ){y ® y) + I(S(y ® y)) 
= I(S 2 )(y ® y) + /(Sy ® <5y) + /(<5y (8) Sy) + /(<5y ® Sy), 

where we have used Lemma 6.1. Expanding 5y in this equation and invoking 
relation (118), we thus obtain 

5z = X~(y ® y) + X^{y v ® y) + X*(y ® y v ) + X v - V (y v ® y v ) 

+ AT^y* ® y) + X*(y^ ® y) + A*(y ® y^) + X^(y ® y ') + z\ 

where z has to be understood again as a remainder. Now, by our standard, 
argument we shall define z° in the nonsmooth case by = — A(J), where J 
is given by 

J = 5[X"{y ® y) + X x '(y v ® y) + X^(y ® y v ) + A x - y (y v ® y v ) 

(121) 

+ X^y^ ® y) + X^(y^ ® y) + X*(y ® y v ) + A^(y ® y ')]. 

In order for this equation to be well defined, we still need to check that 
J belongs to C^ +3k (B- p ), which is in the domain of A since 7 + 3k > 1. Let 
us then compute J: owing to (116), we have 

5[X-{y ® y)] = [6X~](y ® y) - X~[5(y ® y)} 

= -X v (Sy ® Sy) - X v (Sy ® Sy) - X v (<5y ® 5y), 

thanks to relation (113) and Lemma 6.3. The other terms can be computed 
along the same lines, and here is a sample of what is obtained: 

6[X*(y" ® y)] = X v (X v y v ® Sy) - X*(Sy" ® Sy) 

- X^(Sy v ® <5y) - X^(5y v ® <5y), 
8[X^(y* ® y)] = VilCy* ® Sy) + A^A^ ® S 2 )(y x ' ® y) 
-X^y^y) 

and 

<5[X V - V (y v ® y v )] = A v (X v y v ® X v y v ) + A^(X v y v ® S 2 y v ) 
+ X x -(S 2 y v ® A v y v ) - X x - v (S 2 y v ® <5y v ) 
- X v - y (<5y v ® S 2 y v ) - A v - V (5y v ® 5y v ). 
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Now, by gathering all the terms we have obtained, we obtain that J = 
Tl=i J k, with 

Ji = X v [S*y <g> (Sy + A v y v + X^y^ + X^y*) 

+ (-Sy + X v y v + X*y* + X^y*) ® Sy 
- <5y v ® <5y v + (X v ® X v )(y v ® y v )], 
J 2 = X^[-Sy~ ®5y~ <fy v <8 Sy + (5 2 X v )(y v ® y v ) - <5y v ® 5y 
+ (A v ® «S 2 ) (y* ® y) + (S ® X v ® S){y* ® y)] 

and 

J 3 = X^f-Sy ® <5y v -5y® Sy~ + (X v ® S 2 )(y v ® y v ) - <5y ® <5y v 

+ (<S 2 ® X v )(y ® y*') + (S <g) X v ® 5)(y y*)], 

J 4 = -X x - y 5(y v ® y v ) - X^y^ ® y) - X*<%^ ® y) 

- X*5(y ® y v ) - X^5(y ® y '). 

Furthermore, notice that, using equation (119) for the increments of y, the 
quantity J\ can be simplified into 

Ji = X v [-,Sy ® y fl - y fl ® Sy - 5y^ ® <5y v + (X v ® X v )(y v ® y v )]. 

We are now left with the cumbersome task which consists in analyzing 
the regularity of all the terms we have produced so far. We shall just focus 
on one particular example, namely X x '(5y v ® Sy), leaving the other ones to 
the patient reader. Invoking again Lemma 6.1, we have 

X*(6y v ® Sy) = X^(Sw ® 5w ® Sy) + X^(Sw ®5w® Sy) 

+ X^(5w®Sw®Sy). 

Among the three terms in the right hand side of this relation, we shall ana- 
lyze the first one, the other ones being similar: recall that X^ G C^ +2k °C(B^, 

B- p ), Sw G C^(B V ) and Sy G C%(B V ). Thus X*(5w®5w® Sy) G C^ +2K ° +ti (B^ p ), 
which is enough regularity to apply the A-map. The other terms in the de- 
composition of J can be treated similarly, which ends the proof of our first 
assertion. 

Our second claim being immediate from the construction of our integral, 
let us say a few words about the last one. Here again, many terms have to 
be estimated, and we shall focus on a representative example, namely the 
term = y v = w ® w. In fact, the quantity Af[w ® to;C* 2 (y?S)] has to be 
estimated, and recall that, according to Lemma 6.1, the following relation 
holds true: 

5(w ® w) = Sw ®6w + 6w<g) Sw + Sw ® 8w. 
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Thus, since w € CJ" 1 (B v ) and St is a bounded operator on B r] for any positive 
5, we obtain 

M[w ® Cf ([0, S}; £ 2 )] < c(l + |w |„ + S".A/>; Cf ([0, S]; H 2 )]) 2 

< c(l + |wo|„ + S»N[y; QxA% S])])\ 

where we used also the decomposition wt = StoWQ + Swto to bound wt in 
terms of wq and 

Kir, < \wo\ v + S?Af[w; C^([0, S];B 2 V )]. 

The other terms defining Sz can be treated along the same lines, which 
proves relation (120). □ 

We can turn now to the main goal of this section, which is to get an 
existence and uniqueness result for equation (107). 

Theorem 6.8. Assume that x allows to define some incremental opera- 
tors X T for any r G T such that d(r) < 4, and that these increments satisfy 
Hypothesis 4, for j,kq,k,t] such that k < kq < 7, 7 + 3k > 1 and rj > 1/4. 
Then there exists a strictly positive Tq = Tq(X t ; d{r) < 4) such that equation 
(107) admits a unique solution y € Qk,x([0,To]). 

Proof. The proof of this result is very similar to those of Theorems 
4.3 and 5.6, and we shall omit the details here. Just notice that inequality 
(120) allows to construct an invariant ball for the map T in Q k ,x([0,Tq]), 
whenever To is small enough. The contraction argument can then be written 
in a standard way. □ 

6.4. The Brownian case. In this section, we investigate the behavior 
of the operators X T defined above, when x = X is an infinite-dimensional 
Brownian motion, defined at Section 3.4. Our aim is of course to show that, 
under certain conditions, X satisfies Hypothesis 4. To this purpose, for the 
remainder of the section, we will mainly consider some applications on the 
space Bri for 77 = 1/4 + e and a small e > 0. Let us also introduce an addi- 
tional notation: for the remainder of the article, we will write A < B for two 
real quantities A and B when A < cB for a universal constant c. 

Proposition 6.9. Let X be an infinite- dimensional Brownian motion 
defined by the covariance structure (45), with Q given by (44) f or v G 
(1/3, 1/2]. Recall also that, for a planar binary tree, we have set |r| = d{r) — 
1. Let X T be the incremental operator given by (109) where the stochastic 
integrals have to be understood in the ltd sense. Then, almost surely, X T G 
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Cg Ais(£>»7 j^r?) /or t = V) ^, >y, v.y, \,, ^, X', / and for any kq satisfying 
0<k <1/4-?? + ^/2. Moreover, X T e ^ +/to(|r| " 1)_1/4 £Hs(Bj (T) ; B- p ) for 
7 = Ko + r ? + / 5 < 1/2- Theorem 6.8 can then be applied in this situation. 



Proof. In all the cases, the line of the proof is the same. We obtain an 
I? estimate on the Hilbert-Schmidt norm of X T , which by Gaussian tools 
can be boosted to an L p bound for any p. Applying Lemma 3.8, the result 
is then easily deduced. 

Admitting for the moment the results of Lemma 6.11 below, let us give 
some details about our method. Since our incremental processes always be- 
long to a finite chaos of the infinite-dimensional Brownian motion X, it is 
easily deduced from Lemma 6.11 that 

E[\\XT\\ P ... }<(t- a )PM(«0+e+(l-l/M)/2) < , t _ s y\r\(K +e) 

for any 0< kq<1/4 — rj + u/2 < rj/2. Moreover, it also holds that: 

ll«£JI<E ll^llll^ll<(i-^ |rl|K0 («-^ |r2|K0 <(i-^) MK0 5 

where t 1 ,t 2 denote the trees appearing in the expansion for SX[ US given at 
Lemma 6.3, for which we have always |t |, |r 2 | > 1 and | -7 -1 1 + |r 2 | = |r| + 1. 
So it is clear that using the extended G-R-R Lemma 3.8, we obtain 

\\xLs\\<(t-s) lTlK0 , 

for any r such that d(r) < 4. Finally, observe that the conditions 7 < 1/2, 
Kq = u/2 — e and 3kq = 7 > 1 force us to choose u > 1/3. □ 



An easy consequence of the last estimations is an existence result for a 
Brownian SPDE in the rough-path sense: 



Theorem 6.10. Let X be an infinite- dimensional Brownian motion on 
[0,T] x [0, 1], defined by the covariance function given by (45) and (44) with 
u > 1/3. Then there exists rj > 1/4, 0<k<7<1/2 satisfying k < kq and 
7 + 3k > 1 such that, for any ip € B v the equation 

Y(0, = V>(£), d t Y(t, £) = AY(t, dt + Y(t, 2 X(dt, d£), 

t€[0,T],£€[0,l], 

with periodic boundary conditions, understood as equation (84), has a unique 
local solution in Q K ,n,ip U P t° a ^ me ^* which depends on the initial condition 
and on the operators X T , |r| < 3. 
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Proof. Like in the proof of Theorem 5.8, the proof amounts to check 
the validity of Hypothesis 4 in the light of Proposition 6.9. □ 

The rest of the paper is dedicated to the L 2 estimations for the operators 
X T . In fact, we will obtain a slightly stronger result than the one we claimed 
at Proposition 6.9. 

Lemma 6.11. For the trees considered at Proposition 6.9, we have the 
following L 2 bounds: 

(122) £[||A7J 2 , „ w J < (t-.s) |T|A " 1/2 , 

where A = 1 — 2r] + v — e for some arbitrary small e > 0. 

Proof. It is conceptually easy to generalize the arguments of Propo- 
sition 5.7 to reduce the problem to an estimation of a mixed sum (over 
eigenvalues of A Q ) and integral over time variables (after contraction of the 
stochastic integrals). This long and tedious task is left to the reader. We 
prefer to give a diagrammatical algorithm which allows to go from the ker- 
nel on L 2 (associated to each operator) to a simple sum estimation. This 
will be detailed in the next two subsections. □ 

Remark 6.12. We can interpret this result by the following heuristic 
considerations. The situation more similar to the finite dimensional theory 
is when rj is slightly larger than 1/4 and v slightly larger than 1/2. In this 
case, A is arbitrarily near to 1 which would give the classical scaling of 
Brownian increments if we could ignore the factor —1/2 apperaring at the 
exponent in the r.h.s. of equation (122). This further loss of regularity is 
due to the need of estimating the Hilbert-Shmidt norm. We conjecture that 
some technique which would allow to estimate directly the operator norms 
of the increments would give a better time regularity which would improve 
the overall theory. Apart from this technical difficulty there is an intrinsic 
departure from the Brownian regularity due to a loss in space regularity 
which must be compensated via a transfer from time to space regularity 
(allowed by the convolutions). This loss in space regularity has two sources: 
one is the nonlinear operation which start to be badly behaved when ry is 
smaller than 1/4, the other is the presence of the noise which degrades the 
spatial regularity of the result when v is smaller than 1/2. Both contributions 
are clearly accounted in the formula for the effective time regularity A. 

Remark 6.13. As we noted also elsewhere in the context of rough-paths 
associated to deterministic PDEs [7-9], in the infinite-dimensional setting 
objects like the semigroup S (and in general the unbounded linear operators 
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appearing in the equations) must be considered at the same level of the 
driving stochastic processes in the sense that the ensemble of these objects 
form a rough-path. In this perspective, the fact that the rough path X T 
which we construct depends both on the Gaussian noise and on the specific 
convolution semigroup S should not be considered more unusual that the 
fact that in the finite-dimensional theory the higher order iterated integrals 
depends on the vector of all irregular components driving the equation. 



6.5. Diagrammatica. We will first show, at a heuristic level and on a 
simple example, how to pass from an incremental operator to a graph for 
the computation of Hilbert-Schmidt norms. 

(1) Case of the operator X v . Consider first the operator X W :B 2 B v . 
Recalling the notation of Section 3.4, an orthonormal basis for B^ is given by 
j > 1}, where e~j = XJ v ej. Furthermore, the particular form (44) we have 
assumed on the covariance function Q implies that x can be decomposed as 



where {/3 p ;pG Z} is a sequence of independent Brownian motions. Hence, 
setting (•, •) for the inner product in L 2 (S), the matrix elements of X s * are 
given by 

[X? s ] i)jk = (ei,X? s (ej ®e k )) = (ef, J Studx u (S us ej)(S us e k )\ 



{Stu&i] dx u (S us ej)(S us e~k)} 
= E^^Ap^ 2 / df3'P{Stu(ii,e p (S us e j )(S us e k )) 

= 12 WV^ 2 J* e- x ^- u) ^ {u - s) - x ^ u - s) d^. 

p;i=p+j+k s 

Thus, the Hilbert-Schmidt norm of X^ s in C(B 2 ;B ri ) can be written as 



From this simple computation, the following rules appear: 

• Some multiple sums (involving terms of the form A") with constraints on 
the indices appear, due to the fact that {ef,i € Z} is the trigonometric 
basis of L 2 (S). 
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• Some contractions in the sums take place, because of the Brownian stochas- 
tic integrals. 

With the above considerations in mind, we can associate to X^ s the following 
graphical representation: 

i p 




where the solid lines represent factors of S. This is a bookkeeping device for 
the relation between the various indices and time parameters. The computa- 
tion of iSfll-X^Hgg C ^ B 2. B )] corresponds to putting side by side two specular 
copies of this graph and connecting the corresponding top and bottom lines 
(to compute the HS norm), while contracting in all the allowed ways the 
dashed lines (to compute the contractions of the stochastic integrals). Do- 
ing so, we obtain the graph 

u 



E \\ X ts\\nS,C{Bl;B rl ) = J du U 

where we use the following convention: solid lines correspond to factors of 
S, time parameters are attached to vertices, crossed solid lines correspond 
to factors of A 2r, S (coming form contraction of output lines), crossed double 
lines correspond to factors of A~ 2r) S (coming form contraction of input lines), 
dashed lines are associated to factors of Q (coming from the Ito contraction 
of the noise). We fix an orientation for each edge of the graph and associate 
an index to each oriented edge. To each vertex corresponds a constraint that 
the sum of indexes of incoming edges minus indexes of outgoing edges should 
be zero. According to these rules, the formula of the mean squared norm is 
then 



E 
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and the reader can check that this is indeed the expression for the mean 
value of the HS norm of X]? s . 
Consider now the expression 

i+j+k+l=0 

We trivially have 

i+j+k+l=0 

Furthermore, setting q = k + I, applying Lemma 4.6, and recalling that we 
have chosen 77 > 1 /4, we can bound A as 

i+j+q=0 

where we used the relation J2k+i=q X~k' 2n X l 2q < Xq 2v ■ Moreover, assuming 
that v < 2rj we can use again Lemma 4.6 to get 

A<Y Af ^e- 2 ^-) <(t- «)- £ \ 2 r v -\ 

i i 

and choosing a = 1/2 + 2r/ — v + e so that the sum is convergent, we obtain 
< (t — u) A ~ 3 / 2 , with A = 1 — 2i] + v — e. So we proved the graphical equation 




(123) u <( t _ n )A-3/2 



and hence, if we suppose that A — 3/2 > —1, that is, A > 1/2, we have 
obtained that 

E[\\Xr s \\ls, m A) ] < [ du (t - u) A ~^ <(t- s) A ^ 2 , 

which is the desired bound for Lemma 6.11. Let us say a few words about 
the condition A > 1/2: if r] = 1/4 + e, then one has A = 1/2 + ^ — e + e. 
This means that the condition A > 1/2 can be met as soon as v > 0, which 
simply rules out the possibility of considering a space-time white noise at 
this stage. 

(2) Case of the operator X^. With the same kind of considerations as for 
AT V , it can be shown that the the matrix elements of the operator X s ** are 
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,»\i.w= E E ^k 2 ^v ,2 k v/2 

p+j+n=i q+k+l=n 



X / e 



\i(t-u) 0f e -\j(u-s) 



X J\-^(u-v) df3 g e -X k (v~s) e -X l (v-s)^ 

Thus, its Hilbert-Schmidt norm in jC(B^;B v ) can be written as 



l X tellHS,£(B3;B^) - E 
i,j,k,l 



Y,] |2 
tsJijfeil > 



and the following graphical representation can be associated to this last 
expression: 

i p 



dffi \ dft 3 

J S 



u ^ q 




Thus, for the Hilbert-Schmidt norm of X£, we obtain the graph 



e [\\ x ^\\hs,c(b^b v )] 




and the corresponding formula 
£, [||^IIhs,£(B3;B„)] 



2,i 



i+j+k+l=0 



x \-v\-2v\-2r) e -2\ i (t-u)-2\ l (u-v)-2\„(v-s)-2\ (v-s)-2\ k (u-s) 
n+m+o=l 
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The strategy to control this expression is now straightforward: bounding 
the exponential and performing the time integrations gives, for two positive 
constants a, b such that a + b < 2, 

E [\\ x u\\ns,c(B3-B v )\ 



i+j+k+l=0 n+m+o=l 



Using the fact that 2r\ > 1/2, we can reduce this to the bound 
E [\\ X ^\\RS,C(B'f l ;Br,^ 



< 



i+j'+i+fc=0 n+m=l 

Assuming moreover that v < 2rj, we get 

^[||^ll^( B 3 A) ]<(t- S ) 2 ^ +6 ) £ A^A-A- 2 "A r ^. 

«+j+z+fc=o 



At this point, choose 6 = 2rj — u so that i/ + 6 = 2?7 and 

Hence, this sum is finite if we choose a = 1/2 + 2rj — v + e, and we get 

Before proceeding to the estimation of the other more complex operators, 
let us make a useful observation. Consider the following subgraph on the left 
of the previous graph: 




After reduction of the crossed double lines (carrying the factors due to A~ v ) 
by an iterated application of Lemma 4.6, we obtain the following expression 
which corresponds to a bound for this graph: 



< / dv\7 v e~ 2Xi{u ~ v) < 



dv X7U ^ b 



(u — v) 



b v i 
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so that choosing b = 2rj — u, we get an estimate of the form < (y — s) A X i 2r> . 
Summarizing these considerations, we have obtained the graphical equation 



(124) 




which we will use multiple times below. 



6.6. More complex graphs. The tools we have introduced so far will allow 
us to treat the two remaining cases we are left with, namely and X vv . 

(1) Case of the operator X***. By using the same kind of arguments as in 
the previous subsection, we obtain a representation of the form 

P 




Now, invoking repeatedly relation (124), we can iteratively reduce the above 
graph to obtain a bound for ^[11X^11^ c ^ B 4.j 3 )] of the form 

it 




ROUGH EVOLUTION EQUATIONS 73 
< f du (t - u)^l\u - sf A <(t- S ) 3A -V2, 

Js 

which is again what is needed for our Lemma 6.11. 

(2) Case of the operator X v v . Using the same conventions as before, the 
operator X vy can be represented as: 



P 




Now our current situation is slightly different from the previous ones, since in 
the triple Brownian integral above, the last two are not iterated. This means 
that we have to handle some sums of the form E[(^2 Y a Z a ) 2 ] for some 
centered Gaussian random variables (Y a ,Z a ) a forming a Gaussian vector. 
The standard way to compute such sums is to write 

Y,E[Y a Z a Y p Zp] 

a,/3 

E[Y a Z a ]E[YpZp] + E[Y a Yp)E[Z a Zp) 

a,/3 

+ E[Y a Zp\E[Y p Z a ]. 



E 



y ] Y a z a I 



By extrapolating this elementary consideration to our situation, this im- 
plies that, in the computation of £^|jX^- v ||^ s y three different kind of 
contractions are involved. Hence, we also obtain three different graphs: 
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Observe that the first of those graphs already appeared in the study of X^. 
We will then focus on the two other ones. 

The analysis of the second graph above can be started by reducing again 
the crossed double lines, which gives a new graph of the form 




However, at this point, we cannot proceed as in the previous cases, with 
a sequence of reduction of subgraphs, in order to prove the convergence. 
Indeed, in the current situation, some irreducible triangular structures like 



ROUGH EVOLUTION EQUATIONS 



75 



the following appear, at the top and bottom of the graph (126): 



(127) 



dw 3 




where the dotted lines stand for the remaining part of the graph, and where 
we put explicit indexes on outgoing lines and on the edges of the subgraph. 
Note that by the rules we have imposed on the graph, the constraint i + 
j + k = (we consider the dotted lines directed inwards) holds true. The 
contribution of the triangular structure is thus given by 



Si 



i+j+k=0 



■V}) — X q -i(u— w) 



and we can bound this last expression by 

dw 



^ S i+3 



+k=0 



(u — w) 1 



E x -i/ x -6/2 x -6/2 
A g+k A q A q-i ' 



for some b £ (0, 1). The latter sum is finite when u + b > 1/2, which means 
that, by choosing 6 = 1/2 — v + e, the triangular structure yields a bound of 
order 



+ j + k=0{U — s 



,1-6 



Summarizing the previous discussion, it is now easily seen that the structure 
(127) behaves like a simple vertex, up to an appropriate factor of (u — s): 



dw 3 




' u ' 



k <(u-s) A J 
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Using this fact, we can reduce our graph (126) to the following simpler 
structure: 



< 



du{u — s 



,2A 




(t-s) 



3A-1/2 



where the last bound has been obtained similarly to (123). 

Finally, let us associate a bound to the third graph in (125). First, notice 
that this third graph can be reduced to 



t ru ru 

du dv dw 




(128) 



Furthermore, this graph contains the irreducible subgraph 



3-1 + 



U rU 

dv / dw 




which corresponds to the expression 



ru ru 

D = dv dw 

J s J s 

x y^expj— Xj(u — v ) — \k(u — w) 

j+k=i q,l 
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- \j-l+q(u -V)~ X k - q+ l(u - W) 

— \j-l(w + v — 2s) — Xk^ q (w + v — 2s)} 

Introducing an additional parameter b and bounding the exponential terms 
as usual gives 



(u — v) b J s (u — w) 1 



j+k=i q 



,1 A j A k\ A q A j-l A k-c, 



Now, using Lemma 4.6, we can bound the sums over q and I in order to 
obtain 



D< 



u dv f u dw 



y 1 



(u - v) b J s (u - w) b .j^ = . \ b +»\^ 

Thus, choosing b = 2rj — v, we end up with 

D<(u-s)^\f\ 
which means that we have obtained the graphical inequality 




< 



(u-s) 



2A 



Plugging this representation in the complete graph (128), we obtain 

u 



< 



du{u — s) 




< 



(t-s) 



3A-1/2 



Going back to Lemma 6.11, we should still treat the case of X T for r = 
^, ■¥ . But these estimates are now mere variations of the previous ones, 
and are left to the reader for sake of conciseness. 



78 



M. GUBINELLI AND S. TINDEL 



REFERENCES 

[1] Adams, R. A. (1975). Sobolev Spaces. Pure and Applied Mathematics 65. Academic 

Press, New York. MR0450957 
[2] Coutin, L. and QlAN, Z. (2002). Stochastic analysis, rough path analysis and 

fractional Brownian motions. Probab. Theory Related Fields 122 108-140. 

MR1883719 

[3] Dalang, R. C. (1999). Extending the martingale measure stochastic integral with 

applications to spatially homogeneous s.p.d.e.'s. Electron. J. Probab. 4 no. 6, 29 

pp. (electronic). MR1684157 
[4] Feyel, D. and de La Pradelle, A. (2006). Curvilinear integrals along enriched 

paths. Electron. J. Probab. 11 no. 34, 860-892 (electronic). MR2261056 
[5] Friz, P. and VlCTOlR, N. Multidimensional dimensional processes seen as rough 

paths. Cambridge Univ. Press. To appear. 
[6] Gubinelli, M. (2004). Controlling rough paths. J. Funct. Anal. 216 86-140. 

MR2091358 

[7] Gubinelli, M. (2006). Rough solutions of the periodic Korteweg-de Vries equation. 
Preprint. 

[8] Gubinelli, M. (2006). Rooted trees for 3D Navier-Stokes equation. Dyn. Partial 

Differ. Equ. 3 161-172. MR2227041 
[9] Gubinelli, M. (2008). Abstract integration, combinatorics of trees and differential 

equations. Preprint. 

[10] Gubinelli, M. (2009). Ramification of rough paths. J. Differential Equations. To 
appear. 

[11] Gubinelli, M. (2009). Rough integrals in higher dimensions. Unpublished 
manuscript. 

[12] Gubinelli, M., Lejay, A. and Tindel, S. (2006). Young integrals and SPDEs. 

Potential Anal. 25 307-326. MR2255351 
[13] Lejay, A. (2003). An introduction to rough paths. In Semmaire de Probabilites 

XXXVII. Lecture Notes in Math. 1832 1-59. Springer, Berlin. MR2053040 
[14] Leon, J. A. and San Marti'n, J. (2007). Linear stochastic differential equations 

driven by a fractional Brownian motion with Hurst parameter less than 1/2. 

Stock. Anal. Appl. 25 105-126. MR2284483 
[15] Lototsky, S. V. and Rozovskii, B. L. (2006). Wiener chaos solutions of linear 

stochastic evolution equations. Ann. Probab. 34 638-662. MR2223954 
[16] Lyons, T. and QlAN, Z. (2002). System Control and Rough Paths. Oxford Univ. 

Press, Oxford. MR2036784 
[17] Lyons, T. J. (1998). Differential equations driven by rough signals. Rev. Mat. 

Iberoamericana 14 215-310. MR1654527 
[18] Maslowski, B. and Nualart, D. (2003). Evolution equations driven by a fractional 

Brownian motion. J. Funct. Anal. 202 277-305. MR1994773 
[19] Nualart, D. (1995). The Malliavin Calculus and Related Topics. Springer, New 

York. MR1344217 

[20] Pazy, A. (1983). Semigroups of Linear Operators and Applications to Partial Dif- 
ferential Equations. Applied Mathematical Sciences 44. Springer, New York. 
MR710486 

[21] Perez-Abreu, V. and Tudor, C. (2002). Multiple stochastic fractional integrals: 
A transfer principle for multiple stochastic fractional integrals. Bol. Soc. Mat. 
Mexicana (3) 8 187-203. MR1952159 



ROUGH EVOLUTION EQUATIONS 



79 



[22] Quer-Sardanyons, L. and Tindel, S. (2007). The 1-d stochastic wave equation 
driven by a fractional Brownian sheet. Stochastic Process. Appl. 117 1448-1472. 
MR2353035 

[23] Tindel, S. (1997). Stochastic parabolic equations with anticipative initial condition. 

Stochastics Stochastics Rep. 62 1-20. MR1489179 
[24] Tindel, S., Tudor, C. A. and Viens, F. (2003). Stochastic evolution equations 

with fractional Brownian motion. Probab. Theory Related Fields 127 186-204. 

MR2013981 

[25] Walsh, J. B. (1986). An introduction to stochastic partial differential equations. In 
Ecole D'ete de Probabilities de Saint-Flour, XIV — 1984- Lecture Notes in Math. 
1180 265-439. Springer, Berlin. MR876085 

CEREMADE Institut Elie Cartan Nancy 

Universite de Paris-Dauphine Universite de Nancy 

Place du Marechal De Lattre De Tassigny B.P. 239 

75775 Paris cedex 16 54506 Vandceuvre-les-Nancy Cedex 

France France 

E-MAIL: massimiliano.gubinclli@ceremadc.dauphine.fr E-MAIL: tindel@iecn.u-nancy.fr 



